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We construct holographic duals to A/" = 1 SQCD with a quartic superpotential and un- 
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C^ \ dimensional submanifolds of an internal space. The flavor degrees of freedom are intro- 

duced by means of D5-branes extended along two-dimensional calibrated surfaces, and act 

cn \ as sources of the different supergravity fields. The backgrounds we get include the back- 

t^^ \ reaction of the flavor branes and generalize the geometries obtained so far to the case in 

which the fundamental matter is massive. The supergravity solutions we find are regular 
everywhere and depend on a radial function which can be determined from the distribution 
of flavor branes used as sources. We also work out the holomorphic structure of the model 

^ ■ and explore some of its observable consequences. 
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1 Introduction 

The AdS/CFT correspondence [Ij is one of the greatest conceptual developments in the 
study of the dynamics of gauge theories of the recent years (see [2] for a review). This 
correspondence, which is in agreement with early ideas of 't Hooft |3] on the large Nc limit 
of QCD, has provided a whole set of analytical tools to study gauge theories in the strongly 
coupled regime. Besides, through the holographic principle of quantum gravity, the corre- 
spondence has established a fascinating and far-reaching connection between gauge theories 
and the physics of black holes. 

Although the final goal would be finding a dual to QCD, a less ambitious and more feasible 
objective is trying to construct gravitational backgrounds dual to minimal supersymmetric 
models. In this paper we will concentrate on the model [1] dual to J\f = 1 super- Yang-Mills 
(SYM) theory, which is based on the geometry obtained by Chamseddine and Volkov in ref. 
[5] . This geometry can be regarded as generated by a set of N^ D5-branes wrapping a compact 
two-cycle of a Calabi-Yau cone. If the size of the cycle is small, the low-energy description 
of the wrapped D5-branes is effectively (3-|-l)-dimensional. However, at larger energies the 
Kaluza-Klein modes on the compact cycle show up and mix with the four-dimensional field 
theory degrees of freedom. Nevertheless, the model of [1] successfully encodes confinement 
and chiral symmetry breaking in a geometric setup (see [6] for various reviews). 

Another step in the process of approaching the holographic gravitational models to phe- 
nomenology is the addition of flavor, i.e. of fields transforming in the fundamental repre- 
sentation of the gauge group (quarks). From the point of view of the gravitational theory, 
adding quarks to a given gauge theory corresponds to incorporating additional branes to the 
setup \7\. These flavor branes should extend along the gauge theory directions and wrap a 
non-compact cycle in the internal manifold, in order to make their worldvolume symmetry a 
global flavor symmetry. If the number Nf of flavor branes is much smaller than the number 
Nc of color branes, one can reasonably neglect the effect of the flavor branes on the geometry 
and treat them as probes. This is the so-called quenched approximation which, on the field 
theory side, amounts to considering the quarks as external non-dynamical objects that do 
not run in the loops. For the model dual to A^ = 1 SYM the flavor branes are also D5-branes 



which wrap a non-compact submanifold of the internal space in such a way that A/" = 1 is 
preserved [8] (see [9] for a review of similar studies in several other models). 

In this paper we are interested in studying unquenched flavor in the holographic dual of 
A/" = 1 SYM. In this case one has to compute the backreaction by solving the equations 
of motion of a system of gravity with brane sources. Generically, these sources modify the 
Einstein equations and the Bianchi identities of some Ramond-Ramond field strengths. We 
will follow the approach initiated in ^Q\, in which one has a large number of flavor brane 
sources which are delocalized and one has to deal with a continuous smeared distribution 
of branes (see [H] for an earlier implementation of this idea in the context of non-critical 
string theory). In this approach the sources do not contain Dirac 5- functions, which greatly 
simplifies the task of solving the equations of motion. On the field theory side this setup 
corresponds to the so-called Veneziano limit, in which both Nc and Nf are large but their ratio 
is kept fixed. In refs. [Ill [13], [H] different aspects of the supergravity duals of A/" = 1 SYM 
with smeared flavor branes were studied, whereas this approach has been also successfully 
applied to other types of backgrounds (see [I5] for a detailed review). 

The A^ = 1 flavors added in refs. [TOl [121 [13] are massless, which amounts to considering 
flavor branes extended along the full range of the holographic coordinate r. The correspond- 
ing supergravity solutions are singular in the IR. This is, actually, a common feature of most 
massless flavored solutions found so far with the smearing technique (see, for example, refs. 
[ini HTl [18] for the D3-D7 systems on the conifold). This curvature singularity can be quali- 
tatively understood as due to the fact that, for massless flavors, all branes pass through the 
origin r = and, therefore, the brane density is highly peaked at r = (an exception to this 
behavior is the solution recently found in [19] for the gravity dual of Chern-Simons-matter 
theories with flavors). 

To remove the IR singularity one can consider massive quarks or, equivalently, a family of 
flavor branes which do not reach the origin (another possibility is to add temperature and 
to hide the singularity behind a horizon, as was done in ref. [SO])- For the D3-D7 system 
these regular solutions for massive flavors were found in refs. [2 n [221123] . As argued in [16j, 
passing from the massless to the massive case in these systems just amounts to substituting 
in the ansatz Nf by NfS{r), where S{r) is a profile function that interpolates between zero 
in the IR and one in the UV. To calculate S{r) one has to perform a microscopic calculation 
of the flavor brane charge density, whose result is not universal since it depends both on 
the characteristics of the unflavored system and on the particular family of flavor brane 
embeddings. 

In this paper we find supergravity backgrounds dual to A/" = 1 SYM theories with un- 
quenched massive quarks. The first step in our analysis will be finding the precise deforma- 
tion of the background which corresponds to the backreaction induced by massive flavors. 
We will show that the compatibility with the Af = 1 supersymmetry implies a certain type 
of deformation which is also parameterized by a profile function S{r). When this function 



S is identically equal to one we recover the results of [TO] for massless quarks. However it 
is important to point out that, in this D5-brane case, the massive quark ansatz cannot be 
recovered by performing the Nf — )■ NfS{r) substitution in the massless ansatz of [TO] . 

From our ansatz we will be able to obtain a consistent system of first-order BPS equations 
which can be partially integrated and reduced to a second-order master equation which is 
the generalization to this massive case of the equation derived in [13] for massless quarks. 
To solve this master equation (and the full BPS system) one needs to know first the profile 
function S{r) which, as mentioned above, is not universal and depends on the family of 
embeddings of the flavor D5-branes. Such families are generated by acting with isometrics 
on a fiducial representative embedding. It turns out that only a particular set of these 
families produces a backreaction which is compatible with our ansatz. For this reason we 
must generalize the results of p] and find new classes of supersymmetric embeddings of 
flavor D5-branes. In order to carry out this analysis we will introduce a convenient set of 
complex coordinates suitable to represent the metric and forms of the 5'f/(3)-structure of 
our geometry. Employing these variables we will be able to find a family of compatible 
embeddings and to compute the corresponding profile function S{r). 

For massive quarks the function S{r) vanishes when r is less than a certain value tq, 
which is related to the mass of the quarks. For r < ro the BPS system coincides with the 
unfiavored one, which corresponds to the fact that the quarks are effectively integrated out 
in this low-energy region. As shown in refs. [THl 113] there exists a one-parameter family 
of solutions of the unfiavored system which are regular at r = 0. Our flavored solutions 
coincide with these in this < r < tq region and, although a potential threshold singularity 
could appear at r = ro, we will show how to engineer brane distributions which give rise to 
geometries that are regular everywhere. 

The rest of this paper is organized as follows. In section |2] we review the basic features 
of the holographic dual to unfiavored A/" = 1 SYM. In section |3] we study the addition of 
massive flavor to the J\f = 1 background and we present our ansatz for the backreaction 
induced by a smeared distribution of flavor branes. In this section we will also present the 
result of the partial integration of the BPS system, as well as the master equation for massive 
flavors. The holomorphic structure of the model is worked out in section |H In section [5] we 
develop a technique to compute the charge distribution function S{r). In this method S{r) 
is obtained by comparing the Wess-Zumino action for the continuous set of branes and that 
of a single representative embedding. By applying this procedure we will discover that not 
all the families of embeddings produce a backreaction compatible with our ansatz. In section 
|6]we find a simple class of compatible embeddings and we compute the corresponding profile 
function. 

The problem of the threshold singularities is analyzed in section [TJ where we show how 
to avoid them and how one can construct regular flavored backgrounds. In section |8] we 
integrate numerically the master equation and we provide numerical solutions for the different 



functions of the ansatz. Some observable consequences of our model are analyzed in section|9l 
In section [10] we summarize our main results and we discuss some further lines of research. 
The paper is completed with several appendices. In appendix |X] we study in detail the 
realization of TV = 1 supersymmetry for our ansatz and we analyze the corresponding BPS 
system. In appendix [B] we write in detail the equations of motion satisfied by our solutions. 
Appendix [C] contains a microscopic calculation of the charge density for some embeddings. 
Finally, in appendix [D] we reconsider the Klebanov-Strassler model with unquenched massive 
flavors and we apply the new techniques developed in the main text to compute the D7-brane 
source distribution. 

2 The holographic dual of A/" = 1 S YM 

In this section we will briefly review the supergravity dual to A/" = 1 S YM found in ref . [1] , 
which is based on the four-dimensional supergravity solution obtained in [5]. This super- 
gravity background is generated by Nc D5-branes that wrap a compact two-cycle inside a 
Calabi-Yau threefold. At low energies this supergravity solution is dual to a four- dimensional 
gauge theory, whereas, at sufficiently high energy, the theory becomes six- dimensional. More- 
over, due to a twisting procedure in the compactification, the background preserves four 
supercharges. The corresponding ten-dimensional metric in Einstein frame is given by: 

ds?o = 9sa'N, et [ ^T^dx? 3 + dr^ + e'^ ( dO^ + sin^ ed(j)^) + j {u' - Af 1 , (2.1) 

where $ is the dilaton. In the remaining of the paper we will use units where Qga' = 1. The 
angles 6 G [0, tt] and G [0, 2tt) parameterize a two-sphere which is fibered by the one-forms 
A^ {i = 1,2,3), which can be regarded as the components of an SU{2) non-abehan gauge 
vector field. Their expressions can be written in terms of a function a{r) and the angles 
{9, (p) as follows: 

A^ = -a{r)de, A^ = a(r) sin ^d0, A^ = -cosOdcf). (2.2) 



The w' 's appearing in eq. (12. ip are the SU{2) left-invariant one-forms, satisfying do; 
angles 9, (p and ip 



I tijkOJ^ A w^, which parameterize a three-sphere and can be represented in terms of three 



uj^ = cos'?/'d6' + sin ?/) sin 6^ d0 , 

o;^ = — sin tp d9 + cos tp sin 9 d(p , 

u^ = dip + cos9d4) . (2.3) 

The three angles 9, cp and ip take values in the range 0<6'<7r, O<0<27r and () <ip < Air. 
For a metric ansatz such as the one written in (12.11) one obtains a supersymmetric solution 



when the functions a{r), h{r) and the dilaton $ are: 
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where $o is the value of the dilaton at r = 0. Near the origin r = the function e'^^ in (12 ■4p 
behaves as e"^^ ~ r^ and the metric is non-singular. The solution of the type JIB supergravity 
includes a Ramond-Ramond three-form F(3) given by: 

F(3) = -^ {oo' - A') A {u^ - A^) A {oo^ - A') + ^ ^ F-^ A (cD" - A'^) , (2.4) 

a 

where F°- is the field strength of the SU{2) gauge field A"-, defined as: 

F'' = dA- + ]^tabcA''AA'. (2.5) 

When the A'^'s are given by (12. 2p . the different components of F"" are: 

F^ = -a'drAde, F^ = a' sin^dr A d0 , F^ = ( 1 - a^ ) sin^d^ A d0 , (2.6) 

where the prime denotes derivative with respect to r. 

One can readily verify that, due to the relation (12. 5p . the three-form F^^,) written in (12.41) 
is closed, i.e. it satisfies the Bianchi identity dF(3) = 0. Moreover, the field strength (12.41) 
satisfies the flux quantization condition corresponding to Nc color D5-branes, namely: 

i^(3) = N, , (2.7) 



^^10 -'s 



where the three-sphere is the one parameterized by the three angles 9, and ip dX & fixed 
value of all the other coordinates and, in order to check (12. 7p . one should take into account 
that, in our units, T5 = l/(27r)^ and 2/t^Q = (2??)'^. 

It was argued in |3] that the background written above is dual to A/" = 1 SYM in four 
dimensions plus some Kaluza-Klein (KK) adjoint matter. The four- dimensional theory is 
obtained by reducing the six-dimensional theory living on the D5-branes with the appropriate 
topological twist. The latter is necessary to realize the M = 1 supersymmetry on the 
curved space [23] . The KK modes in the four-dimensional theory have masses of the order 
1/a/(7s o.' Nc- Since this mass is of the order of the strong coupling scale, the dynamics of 



the KK modes cannot be decoupled from the dynamics of confinement. A proposal for a 
concrete lagrangian of the vector A^ = 1 multiplet and the different KK modes has been 
written in [13] (see also |25]). Schematically, this lagrangian has the form : 

L = Tt[-^F'^, - tX^D.X + L(<Dfc, Wk, W)] , (2.8) 

where $fc and Wk represent the infinite number of massive chiral and vector multiplets and 
W denotes the curvature of the massless M = 1 vector multiplet V = (A, A^). 

Let us finish this section by recalling that there exists another solution of type IIB super- 
gravity which is directly related to the one written above. In this solution the metric and 
the RR three-form are also given by the ansatz (I2.ip - (I2.6I) but, in this case, the function 
a{r) vanishes, e'^^ = r and e^*~^*" = j^. Actually, these functions are just the UV limit 
(r —7- oo) of the ones written in (12.41) . On the other hand, at r = this new background 
has a (bad) singularity that is solved by the turning on of the function a{r) in the solution 
(12. 4p . which makes the A"' a non-abelian one- form connection with components along the 
three SU{2) directions. This way of resolving the singularity is related, on the field theory 
side, with the phenomena of confinement and R-symmetry breaking of A^ = 1 SYM. 

3 Addition of massive flavors 

Let us now introduce flavors by means of pairs of chiral multiplets Q and Q transforming in 
the fundamental and antifundamental representations of both the gauge group SU{Nc) and 
the flavor group SU{Nf). The lagrangian for the {Q,Q) fields is given by the usual kinetic 
terms and the Yukawa interaction between the quarks and the KK modes, which can be 
schematically written as: 

Lq^q= f d^e {q^ e-"" q + qU"" q) + fd^eq^kQ. (s.i) 

In the effective low-energy theory obtained by integrating out the massive modes, the Yukawa 
coupling between (Q, Q) and the $fe gives rise to a quartic term for the quark fields (see 
[T0l[l2l[T3] for details). 

On the gravity side the addition of flavors can be performed by means of flavor branes, 
which add an open string sector to the unflavored closed string background. For the A/" = 1 
geometry of section [2] the flavor branes are D5-branes extended along a non-compact cycle 
of the Calabi-Yau threefold p]. If the branes reach the origin r = of the geometry, the 
corresponding flavor fields are massless. If, on the contrary, the D5's do not reach r = 0, the 
quark fields are massive (the minimal value of r attained by the brane is related to the mass 
of the quark fields) . 

In this paper we are interested in getting a holographic dual of the A/" = 1 model with 
unquenched matter, in which the dynamics of fundamentals is encoded in the background. 
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To achieve this goal we must go beyond the probe approximation and find a solution of the 
equations of motion derived from an action of the type: 

S = SjjB + Sbranes j (3-2) 

where Sub is the action of ten-dimensional type JIB supergravity and Sbranes denotes the 
sum of the Dirac-Born-Infeld (DEI) and Wess-Zumino (WZ) actions for the fiavor branes. 
Generically, the branes act as sources for the different supergravity fields. In particular, the 
WZ term of Sbranes IS a source term for the RR fields which induces a violation of the Bianchi 
identity of the corresponding RR field strength. In our case, the WZ term of the action of 
a set of D5-branes is: 

Nf 

Swz=n E /,. **(C'(6)), (3.3) 

where C(6) is the RR six- form potential and z*(C(6)) denotes its pullback to the D5-brane 
worldvolume. Let us rewrite (13. 3 p as a ten-dimensional integral, in terms of a charge distri- 
bution four-form Q: 



Swz = n / C(6) A n . (3.4) 

JMio 

The term (13. 4 p induces a violation of the Bianchi identity of -F(3). In order to determine it, 
let us write supergravity plus branes action ( 13. 2p in terms of the RR seven-form F(7) and its 
six- form potential C(6). This action contains a contribution of the form: 

- i ^ / e"* F(7) A *F(7) + Ts / ^6) A n . (3.5) 

^'^lO ^ JMio JMio 

The equation of motion of C(6) derived from (13. 5p gives rise to the Maxwell equation for F(7) 
with Q playing the role of a source, which is just: 

d(e-* *F(7)) = -2KlQT,n. (3.6) 

Taking into account that, in our units, 2k1qT5 = (27r)^ and that F(3) = — e^* * F(7), we get 
that (13. 6p is equivalent to the following violation of Bianchi identity of F(3): 

dF(3) = Air^n . (3.7) 

The four-form Q is just the RR charge distribution due to the presence of the D5-branes. 
Clearly, Q is non-zero on the location of the sources. In a localized setup, in which the 
Nf branes are on top of each other, Q will contain Dirac 5-functions and finding the corre- 
sponding backreacted geometry is technically a very complicated task. For this reason we 
will separate the Nf branes and we will distribute them homogeneously along the internal 
manifold in such a way that, in the limit in which Nf is large, they can be described by a 
continuous charge distribution fl. 



As we will detail below, the continuous set of flavor branes that we will use in our construc- 
tion can be generated by acting with the isometrics of the background on a representative 
fiducial embedding and, therefore, all the branes of the continuous set are physically equiv- 
alent. Actually, we will not choose an arbitrary distribution of branes. First of all, we will 
require that all branes are mutually supersymmetric (and thus they will not exert force on 
each other) and that they preserve the same supercharges as the unfiavored background. 
Moreover, we will also require that the deformation induced on the metric is mild enough, 
in such a way that it reduces to squashing the unfiavored metric (12. ip by means of squash- 
ing functions that depend only on the radial coordinate r. One can prove that the most 
general squashing of this type compatible with the A^ = 1 supersymmetry of the unfiavored 
background is the one in which the size of one of the fibered directions in the metric (12.11) is 
different from the other two. Accordingly, we will adopt the following ansatz for the Einstein 
frame metric of the flavored theory: 

ds^ = e2^('^) [da;?^3 + e^'^^^Mr^ + e^'^^'Xde^ + sin^ 6 d^^) + 

(3 8) 

p2g(r) 2fc(r) ^ ' ^ 

+ —-{{u^ + a{r)d9f + {u^ - a{r) sine dcPf) + (w^ + cos^d^)^ . 

Notice that the ansatz (13. 8p is exactly the same as the one considered in [TU] for the case of 
massless flavors. 

Let us next consider the deformation of the RR three-form F(3). Clearly, due to the 
modified Bianchi identity (13. 7p that must be satisfied in the flavored case, F(3) cannot have 
the same form as in (12. 4p . Actually, we will slightly modify (12. 4p and we will adopt the 
following ansatz for the RR three- form F(3): 



^(3) 



~ {oo^-B')A{u;^-B^)A{u^-B'') + ^ ^ (F'^ + T ) A (cu'^ - 5") , (3.9) 



where 5" is an SU{2) one-form gauge connection and F"" is its two-form field strength, 
defined as in (12. 5p . namely: 

F" = dB"" + - e"''^ B^ AB^ . (3.10) 

In (13. 9p . the /" are two- forms that parameterize the violation of the Bianchi identity and 
thus the flavor deformation of the RR three-form. Indeed, when f°' = the three-form F(3) 
is closed by construction, due to the relation (13.10p between F" and S". We will take, as in 
[To] , the following ansatz for 5°: 

B^ = -b{r)de , 52 = b{r) sin^d(/. , B^ = - cos 6 d(j) , (3.11) 

where b{r) is different from the fibering function a{r) of the metric (they are equal in the 
background of [1]). By applying the definition (I3.10p . we get that the different components 
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of the two-form field strength F"" are: 

F^ = -6' dr A d^ , F^ = h' sin ^ dr A d0 



;i-6^) sin^d^Ad^ . (3.12) 



We will adopt for the flavor deformation two-forms /" an ansatz that parallels F", namely: 

/I = _Li(r)drAd^ , /^ = Li(r) sin^dr A d0 , /^ = L2(r) sin^d^ A d0 , 

(3.13) 
where Li and L2 are two functions of the radial variable to be determined. Actually, after a 
detailed study of the realization of supersymmetry for the metric ansatz (13.81) one can show 
that (I3.13P gives rise to the most general form of -F(3) . By computing the exterior derivative 
of (13.91) and applying (13. 7p . one gets the following expression of the smearing form Vt: 



n 



iVc 



+ 



167r2 



167r2 



sin 6* d6' A d0 A 

dr A d6' A w^ A w^ -h 



LoCJ^ A w^ 



L2 dr A Co^ 



+ 



A 



[ sin duj^ MJj^ + cos 6* d6' A u;^ 



(3.14) 



One can now study the realization oi M = 1 supersymmetry in type JIB supergravity for a 
background with metric and RR three-form given by the ansatz written in (13. 8p and (13. 9p . 
This analysis is performed in detail in appendix |A] and leads to a system of first-order BPS 
equations for the different functions of the ansatz. Combining these equations, a partial 
integration is possible. Let us summarize in this section the results of this study of the BPS 
equations. First of all, one can verify that the functions Li and L2 parameterizing /" and Vt 
are not independent. Actually, from the BPS system one can prove that Li can be written 
in terms of the derivative of L2 as follows: 



Li 



L' 



2 cosh(2r) 



Therefore, if we define the function S{r) as: 

Nf S{r) = 



-N.L^ir 



(3.15) 



(3.16) 



then, the two-forms /" of (13.130 become: 

Nf S\r) 



f 



2Nc cosh(2r) 

N 



dr Ad9 



f 



Nf S'(r) 



2Nc cosh(2r) 



sin ^ dr A 



f = -— ^ S{r) sin^d^A 



(3.17) 
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whereas the smearing form Q can be written in terms of S as: 



n 



167r2 



sin OdOAdcbA 



Su'^ Au^ - S'dr Aoj 



~,3 



+ 



H dr A 

327r2 cosh(2r) 



d9 A OJ^ A u-" + d(h A { sin 9 u' A uj-" + cos 9 d9 A u 



~,1 A ,~,3 



,~,2 



Moreover, the function 6 parameterizing the one-forms i?" can be written as: 

^ 2r + VJr) 
^ ' sinh(2r) ' 

where ri{r) is defined as the following integral involving S: 



(3.18) 
(3.19) 



r]{r) 



Nf 



2N, 



tanh(2r) S{r) + 2 dp tanh2(2p) S{p) 



(3.20) 



It follows from these results that the RR three-form F(3) in (13. 9p is determined in terms of 
a unique function S{r). Notice that the case of massless flavors studied in ^0\ is recovered 
by taking S* = 1 in our formulas. Indeed, in this case only the first term on the right-hand 
side of (I3.18P is non-zero and the charge density distribution Q is independent of the radial 
variable. Moreover, by computing the integral in (I3.20p one can show that our ansatz for 
F(3) is reduced to the one adopted in [TO] . 

In the case of massive flavors one expects the charge distribution to depend non-trivially 
on the radial coordinate and, actually, to vanish for values of r smaller than a certain scale 
related to the mass of the quarks. In our approach this non-trivial structure is encoded in 
the dependence of the function S on the radial variable. Notice also that S should approach 
the massless value S' = 1 as r — > oo since the quarks are effectively massless in the deep 
UV. The way in which the profile function S{r) interpolates between the IR and UV values 
depends on the particular set of D5-branes that constitutes our delocalized source and should 
be obtained by means of a microscopic calculation of the charge density (see below). 

Another interesting observation is that, contrary to the backgrounds with massive flavors 
studied in [211 1221 [23], in this case passing from the massless to the massive case is not 
equivalent to substituting Nf by Nf S{r) in the massless ansatz. Indeed, it is immediate to 
check that making this substitution only the first line in (I3.18P is generated, while the last 
two components of Q (which are essential for the consistency of the approach) are missing. 
Notice that these last two terms in Q are precisely those in (I3.14p which are proportional to 
the function Li which, according to (I3.15p . always vanishes when r — )■ oo. This means that, 
in the UV, the two-forms /" that implement the flavor deformation of F(3) are non-vanishing 
only along the third SU{2) direction, while the other two components are excited when we 
move towards the IR. Interestingly, this structure is reminiscent of the way in which the 
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singularity of the particular solution of the unflavored theory reviewed at the end of section 
[2] is resolved in the full solution (12 .ip - ( 12^61) . namely by turning on the function a{r) and 
making the two-form F"- the field strength of a non-abelian magnetic monopole. 

Actually, as shown in appendix [XJ it turns out that one can also integrate partially the 
BPS system for the functions of the metric in terms of S{r). First of all, the function / is 
related to the dilaton $ as: 

/ = J • (3.21) 

Moreover, the dilaton can be related to the other functions h, g and k as: 

h+g+k 

e-'" = 2e-2*o ^— — , (3.22) 

smh(2rj 

where $o is a constant. In order to represent the remaining functions of the metric (13. 8 p let 
us define, following [T3], the functions P{r) and Q{r) in terms of a and g as: 

Q = (acosh(2r) - l) e^^ ^ P = ae^^ sinh(2r) . (3.23) 

The inverse of this relation is: 

e^^ = P coth(2r) - Q , a = „ u^ ^^ n ■ x.i^ \ ' ^^'^^^ 

P cosh{2r) — C^smh(2r) 

It is demonstrated in appendix |X] that, from the BPS system, one can express h and k in 
terms P, Q and S, namely: 

"" " 4 Pcoth(2r)-g ' "" ~ 2 • ^^-^^^ 

It follows from eqs. (I3.2ip - (l3.25p that the dilaton and the functions of the metric are deter- 
mined in terms of P, Q and S. Actually, the function Q can be integrated in terms of the 
profile S{r), namely: 



Q = coth(2r) 



2i V, - NfSj p) 
coth^(2p) 



^P — ., 2/^ X ^ 9o 



(3.26) 



where go is a constant of integration. Moreover, as in pL3j , one can find a master equation: 



One can first notice that in the case S" = 1 eq. (I3.27P reduces to the equation found in 
Otherwise, knowing the function S (from a microscopic description of the smearing), one 
can get Q from (I3.26P and solve the second-order master equation (I3.27p for P. As argued 
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above, each solution of this equation will give a complete solution of the problem. Moreover, 
in appendix |B] we have explicitly written the equations of motion derived from the type IIB 
supergravity plus sources action. One can check that any solution of the BPS system also 
solves the second-order equations of motion written in appendix [Bl 

Finding an analytic solution of this master equation is probably not possible, but we will 
be able to find numerical solutions, and their asymptotics. In order to achieve this goal 
we will have first to identify a family of supersymmetric embeddings whose backreaction 
on the background is compatible with our ansatz and then we must be able to compute 
the corresponding profile function S{r). In the next section we will start to develop the 
machinery necessary to carry out this calculation. 

4 Holomorphic structure 

As stated at the end of section[3], to find the profile function S{r) we must analyze the families 
of supersymmetric embeddings of the flavor D5-branes. This problem was addressed in ref. 
[8] by looking at the realization of kappa symmetry for probe D5-branes in the unflavored 
background described in section O The analysis of |[8j was performed in terms of the angular 
coordinates of the metric (12. ip and some particularly interesting embeddings were found. For 
our present purposes we clearly need a more systematic approach, which could allow us to 
study different families of embeddings and to determine whether or not their backreaction 
is consistent with our ansatz f l3.8p - fl3.13p . As the internal manifold of our background is 
complex, it is quite natural to work in a system of holomorphic coordinates. The purpose of 
this section is to define these complex coordinates and to uncover the holomorphic structure 
of our background. 

Let us begin by introducing a set of four complex variables Zi (i = 1 . . . 4) parameterizing 
a deformed conifold, i.e. satisfying the following quadratic equation: 

ZiZ2-Z3Z4 = l. (4.1) 

We will also introduce a radial variable r, related to the Zi as: 

4 

"^\zi\^ = 2cosh(2r) . (4.2) 

i=l 

In order to find a useful parameterization of the Zi^s, let us arrange them as the following 
2x2 complex matrix Z: 

Z3 Z2 



z z ■ (4-3) 

-Zi —Zi 

Then, the defining equations (14. ip and (14. 2 p can be written in matrix form as: 

det (^) = 1 , Tr {Z Z^) = 2 cosh(2r) . (4.4) 
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It is immediate to verify that the matrix 

Z. = ( .l-r ) (^-S) 

is a particular solution of f l4.4p . The general solution of this equation can be found by 
realizing that the equations in (14 ■4p exhibit the following SU{2)l x SU{2)^ symmetry: 

Z^LZR\ LeSU{2)L, ReSU{2)R. (4.6) 

A generic point in the conifold can be obtained by acting with isometrics on the point (14. 5p . 
Actually, if we parameterize the SU{2) matrices above in terms of Euler angles as: 

a = cos 2^2 , 

= sm I e 2 , 
k = cos 2^2 , 

R={ ; j: ) (4-7) 



a 


-b 


b 


a 


k 


-J 


I 


k 



I = — sm ^e 2 



(4. 



then, the four complex variables zi, Z2,Z3, Zi^ that solve (14 ■4p are given by: 

zi = —e 2 '-'P^'P' e^ 2 sm - sm e 2 cos - cos ■ 

'22 2 

i((h+i) I r+i^ ^ ^ -r—i^ .0.9 

Zi = e^^^^^' e^ 2 cos - cos e 2 sm - sm ■ 

' 2 2 2 2 

_ ^'^ U _ _V' t/ u 

Zn = e2^'^~'^' I e''^*2- cos - sin - + e '' *^ sin - cos 
^ ' 2 2 2 2 

Z4 = —e~2'^'f~f> I e''"^*2 sin - cos - + e~''~*^ cos - sin - , , 
' 2 2 22/' 

where ip = ipi + ip2- We will show below that these holomorphic coordinates are very 
convenient to analyze the supersymmetric embeddings in our flavored backgrounds. It is 
also useful to introduce a new set of complex variables Wj, related to the Zi by means of the 
following linear combinations: 

Zi + Z2 Z\ — Z2 Z^ — Z4 Z^ + Zi 

^1 = ^—, ^2 = ^—, ^3 = ^—, ^4 = ^—. (4.9) 
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These variables satisfy: 

(wi)' + {W2f + (ws)' + {w^f = 1 , (4.10) 

and there is an obvious SO (4) invariance that is obtained by rotating the Wj's. The so-called 
S'0(4)-invariant (l,l)-fornis are defined as [26] : 

T]i = 6'^dwiAdwj , ri2= {6'^Widwj) A[6^''Wkdwi) , rf'j, = e'^^^WiWjdwkAdwi . (4.11) 

In terms of the radial and angular coordinates these forms are given by: 

?7i = — i ( cosh(2r) dr A (Cj^ + cos6'd0) sinh(2r) [sin^d^ A d0 + sin^d6' A d 

r]2 = i sinh (2r) dr A (cD^ + cos 9 d0) , 

r]3 = -i (- sinh(4r) (^sin 61 d6' A d0 - sin 9d6 A d(pj sinh(2r) [d9 Au^ + sin 9d(j) Au'^)\ . 

(4.12) 
The fundamental two-form J of the 5'[/(3)-structure can be written in terms of the rji forms, 
which will be very useful in what follows. In order to find the corresponding expression, let 
us notice that J has been written in (1A.4P in terms of the angle a that rotates the projections 
of the Killing spinors. By using the value of a written in (1A.15I) and the relations (1A.16I) . 
one easily proves that J can be written as: 

e"2 J = — dr A (u^ + cos9d(j)) H ^VA (d9 Au^ + sm9 d(f) Au^) - 

2 ^ ^4 smh(2r) ^ ^ 

e^^ /a cosh(4r) — cosh(2r) . ^ ,^ , , cosh(2r) — a . ~ ~ 
' ^ ^ ^ • sm 6* d6' A d0 H . \ / ^ — sm 6* d6' A 



sinh(2r) sinh(2r) 

(4.13) 
In terms of the r^j's, one can rewrite J as: 



fl— .„ _ . .2, f^^ ^ ^Qsh(2r) ^^^ ^ ,^ acosh(2r)-l ^ ^^_^^^ 



sinh (2r) V sinh (2r) / sinh (2r) 



e 2 J = — 
2i 

Let us now check that the complex variables Zi defined in (14.81) are good holomorphic co- 
ordinates for the internal manifold. Indeed, since our six-dimensional internal manifold is 
a complex manifold, we can write its metric in terms of the (1, l)-form J. Actually, if one 
writes J as: 

J=^h^^dz'' Adz^ , (4.15) 

which is allowed thanks to the fact that J is a (1, l)-form, then one can prove that the metric 
of the internal space can be written as: 

dsl = -/i„^(d2" O dz^ + dz^ O dz") , (4.16) 
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where we have spht the ten-dimensional metric (13. 8p as ds^ = e^/^dx^g + dsg. The h^^ 
coefficients appearing in fl4.15p and fl4.16p can be read from f l4.14p by using the relation 
between the r^j's and the Zi coordinates (see (14. lip and (I4.9P ). Moreover, by using again 
flA.lSp and flA.16p . one can write the three-form flhoi of flA.4p as: 

^hoi = - . ^,^ e^'^+s+'^+^-dz, A dz2 A dz3 . (4.17) 

smh(2r) Zs 

Furthermore, taking into account (I3.22p . we can write Qhoi as: 

Z Zz 

which shows that ^hoi is, indeed, a holomorphic (3, 0)-form for the complex structure corre- 
sponding to the coordinates (14. 8p . 

The RR six-form potential C(6), defined as -F(7) = — e* * F(3) = dC(6), can also be written 
in terms of the rii one-forms. In fact, it follows from the first of the BPS conditions in (JA.SP 
that C(6) can be written in terms of the form J as: 

qg) = e^ d^x A J , (4.19) 

where d^x = dx" A dx^ A dx^ A dx^. Obviously, since J can be written in terms of the 
z^ variables, the six-form C(6) can also be written as a (1, 1) form in the internal space. 
Notice that C(6) is related to the calibration form of a D5-brane, whose pullback to the 
worldvolume determines if the embedding is supersymmetric or not. Having C(6) written in 
complex coordinates is very convenient from the technical point of view since it will allow 
us to analyze the different supersymmetric embeddings by employing the full machinery of 
the complex variables. 

Another relevant quantity that should be invariant under the 5*0(4) isometry is the smear- 
ing form fi in (I3.18p . since it is giving us the charge distribution of the system. It is a 
(2,2)-form which can be cast in terms of (l,l)-forms as follows: 



2NfS f 2cosh2r \ N j S' ( \_ 

sinh^2r \ sinh^ 2r / sinh^ 2r \ cosh2r 



levr^l] = -—rnr ^^^W^ —ttt^^i ) + — TIHT- ^2 a ( r^i - _^^^J ]z ) . (4.20) 



4.1 Supersymmetric embeddings 

It is now straightforward to show that any embedding defined with holomorphic functions of 
the complex coordinates is supersymmetric. Let us study the case of an embedding extended 
in the Minkowski directions, and defined in the internal space in the following way: 

Z2 = F{zi) , zs = Giz,) , Z2 = F{z,) , zs = G(^i) , (4.21) 
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where, for definiteness, we have chosen zi and zi as worldvolume coordinates in the internal 
space. Recall that z^ = z^'^{zi Z2 — !)■ The calibration form /C for a D5-brane in Einstein 
frame is given by: 

/C = e*d^x A J = e-t C(6) . (4.22) 

By using f l4.15p one can easily get the puUback of this calibration form on the worldvolume 
of the embedding, namely: 

2*(/C) =ie*ird^xAd^iAd% , (4.23) 

where we have defined the function K as: 

K = ^{hu + F'h^-2 + G'h,-3 + F'hri + F'F'hr2 + F'G'h23 + G'hsi + G'F'h,-2 + G'G'hs-3) . (4.24) 

Now, we look at the induced metric ds^ on the worldvolume of the embedding. We get from 
gH: 

dsl = e*/Mx?3 + 2ird^id^i . (4.25) 

Therefore, det^f = e'^'^K^, and one has 



V-det^ d^x A dzi Adzi=i e^K d^x A dz^ A dzi = i* (/C) . (4.26) 

This means that the embedding is supersymmetric, proving explicitly that all holomorphic 
embeddings are supersymmetric. 

5 Charge distributions 

The supersymmetric D5-brane embeddings we are looking at are characterized by two alge- 
braic equations of the type: 

Fi{z,) = , F2{z,) = , (5.1) 

which define a non-compact two-cycle C2 in the internal six-dimensional manifold. As argued 
above, the preservation of supersymmetry is ensured if the two functions in (15. ip are holo- 
morphic. However, in the brane setup we are considering we will not deal with a particular 
embedding of the flavor D5-branes but, instead, with a family of equivalent embeddings. 
This family can be generated from a particular representative of the form (15.10 by acting 
with the SU{2)l x SU{2)r isometries of the conifold. Let us recall how these symmetries 
act on the holomorphic coordinates. Under SU{2)l the holomorphic coordinates transform 
as Zi — )■ £•, where 



'l! 



Z3 ^2 \ ^ / az3 + (3zi az2 + (iz^ , , ^. 

-z\ —Z4 J y —az\ -|- fizz —azi + (5z2 ' ' 
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with |ap+|/3p = l. Similarly, the S'f/(2)^ transformation is: 

^3 Z2 \ f 723 - ^Z2 6z^ + 72:2 



zi —Z4 J y — 7-21 + Sz4 —6zi — 72:4 ' ' 

where the complex constants 7 and 6 satisfy the condition I7P + \6\'^ = 1. We want now to 
determine the charge distribution four-form Q (parameterized by the profile function S{r)) 
for a given family of embeddings. Actually, we will employ a procedure which does not 
require performing the detailed analysis of the whole family and that allows to extract the 
function S{r) by studying one single particular embedding belonging to the family [27]. This 
method is based on the comparison between the action for the whole set of Nf flavor branes 
and the one corresponding to a representative embedding. We can choose to compare either 
the DBI or the WZ part of the actions, since supersymmetry guarantees that they are the 
same. The WZ term of the action of the full set of D5-branes is given by the following 
ten-dimensional integral: 

J Mw 

whereas the action of one of the embeddings is just: 

S7t = T, I f (C(6)) , (5.5) 

JMa 

with M.Q being the worldvolume of the representative embedding chosen and i* (C(6)) denotes 
the pullback of C(6) to A^q- Since all the embeddings of the family are related by isometrics, 
they are equivalent and their actions should be the same. Thus, we should have: 

gsrneared ^ j^^ gsingle _ ^^g^ 

The left-hand side of (15. 6p can be obtained by plugging the expressions of ^2 and C(6) written 
in (13.181) and (14.191) respectively. After integrating over the angular coordinates, one gets a 
remarkably simple expression, namely: 



Sheared ^ 27r NfT^ / d'x dr e'* (e^''S + ]-e^^ tanh(2r)5' J . 



(5.7) 



The non-compact two-cycle C2 that the D5-branes wrap can be parameterized by the radial 
coordinate r and an angular variable. After integrating over the latter, the WZ action (15.51) 
can be represented as: 

S'^T = 271 T5 fd^xdre'^Sir) , (5.8) 



where the function S{r) is related to the integral of the pullback of J along the two-cycle by 
means of the expression: 

/ t*{J) = 2n I dr et S{r) . (5.9) 



Notice that the dependence on the dilaton of the right-hand side of eqs. (15 .Sp and (15.91) is 
consistent with the one displayed in (I4.14p and (I4.19P . By plugging (15 .70 and (15. 8p into (15. 6p 
we arrive at the following relation between the profile S{r) and the function S{r): 

e^'^S + ^ e^3 tanh(2r) S' = S{r) . (5.10) 

The function S appearing on the right-hand side of (l5.1Up depends both on the embedding 
and on the different functions of our ansatz. In the case in which S depends only on the 
functions k and g and this dependence is the same as on the left-hand side of (I5.10p . it is 
possible to obtain the profile function S from (15.1 OP . However, this is a highly non-trivial 
condition which most families of embeddings do not satisfy. To illustrate this fact let us 
consider the families of massive embeddings obtained by acting with the SU{2)l x SU{2)fi 
isometrics on the two non-compact two-cycles found in [8]. The first of this two-cycles 
is the so-called unit-winding embedding (see section 6.1 in [8]), which has the following 
representation in terms of the real coordinates of the metric (13. Sp : 

sinh T ~ "^ 

sinhr = -^ , 9 = 9, = 0, ip = ii , (5-11) 

sine* 

where Vg is a constant. In terms of the complex coordinates (14. 8p . one can easily show that 
(15. lip is a particular solution of the following two holomorphic equations: 

-2^1 ^2 = COsh^Tg , Z^ + Z4 = . (5-12) 

By using (15. lip and (I4.13P it is straightforward to compute the pullback of J and to obtain 
S{r). One gets: 



sinh^ r„ ( ^i. o„ sinh^ r 



S{r) = Wl - '' e^'^ + e^3 a , , coshV . (5.13) 

sinh r \ sinh r — sinh r^ 

Notice that the right-hand side of (I5.13P contains the function a{r), which is not present on 
the left-hand side of (I5.10p . Therefore, the determination of the profile S is not possible in 
this case. Similarly, one can consider the so-called zero-winding embeddings of section of 6.2 
of [8]. In this case the cycle is characterized by the equations: 

• V, ^o ^ sinh(2rg) . ~ cos^ ~ ~ 

smh(2r) = ——- , sm9 = — — — , = 0o , V' = tt , (5.14) 

smc^ cosh (zrg) 

which solve the following system of two complex equations: 

Z1Z2 = - , zi - e-^'^" e"^"^" Z4 = . (5.15) 
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In f l5.14p and flS.lSp 0o is a constant. By computing S{r) in this case one gets: 




Sir) = ,/l - '^^^^ e- + e^^ ,/"^'^'^\^ ^ f 2a cosh(2r) - 

*' sinh^(2r) ^ sinh2(2r) - sinh^(2rg) V ^ 

(5.16) 

which, as in fl5.13p . contains the function a{r) and, as a consequence, it is not of the form 

displayed in f IS.lOp . 

Our interpretation of the fact that S is not of the form f IS.lOp for the embeddings (15. lip and 
(I5.14P is that their backreaction is not compatible with our ansatz. Notice that our F(3) in 
(13. 9p . as well as the charge- density four- form fl in (IS.lSp . are dictated by the S'f/(3)-structure 
of the A^ = 1 supersymmetry and they are highly asymmetric with respect to the exchange 
{9,(f)) 4-> {0,(j)). This interpretation is supported by an independent microscopic calculation 
of fl, which we present in appendix O Indeed, we show in this appendix that a simple 
embedding whose S is not of the form (IS.lOp gives rise to a charge-density Q which does not 
fit into our ansatz. Thus, in order to proceed further with our formalism we have to find 
a concrete example of compatible embeddings and we have to determine the corresponding 
charge profile. 

Fortunately, we have been able to find a simple family of embeddings for which S{r) 
depends on the functions k and g in the same way as the left-hand side of (IS.lOp and, as 
a consequence, one can directly read the profile function S{r) for this configuration. The 
profile S{r) obtained in this way can be used to get Q{r) from (I3.26p . and as an input for 
solving the master equation (I3.27p . In the next section we will present these embeddings 
and we will determine the corresponding profile. 

6 A simple class of embeddings 

As shown above, the massive embeddings found in [8] do not seem to produce a backreaction 
compatible with our ansatz and with its underlying S'[/(3)-structure. In principle, we should 
consider the logical possibility that such a compatible set of embeddings does not exist. In 
this section we will discard this possibility by finding a family of embeddings for which (IS.lOp 
can be solved and a simple expression for the profile function S can be found. We confirm 
this fact in appendix [C] by means of an explicit microscopic calculation in the UV region 
of large r of the charge density four-form Q. By considering the full distribution of fiavor 
branes we will indeed show that, for the embeddings discussed in this section, the density 
Q is of the form displayed in (I3.18P and we will find an expression of S{r) which is just the 
limit of the one found by solving (IS.lOp for large r. 

In terms of the holomorphic coordinates (14. 8 p the simplest embeddings one can think of 
are those characterized by two linear relations of the Zj's. Many of these embeddings are 
related by the action of the SU{2)l x SU{2)ji symmetry and they belong to the same set. 
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Instead of considering the full set we will only deal with a particular representative. By using 
the machinery developed above it is rather easy to consider systematically the different linear 
embeddings, to compute the pullback of the fundamental two-form J and to verify if the 
function S{r) depends on the functions of the ansatz as the left-hand side of ( IS.lOp . Most 
of these linear embeddings do not give rise to a compatible charge density. Let us now work 
out an example for which everything works fine. The representative embedding we want to 
focus on can be written in terms of the holomorphic coordinates (14. 8 p as the following two 
linear equations: 

z:i = Azi , z^ = Bz2 , (6.1) 

where A and B are two complex constants. We can parameterize the two-surface defined by 
(16. ip in terms of, for example, zi. 

11 . B 1 , , 

Zo = : , z-i = Azi Za = ; . (6.2) 

1-ABzi' ^ 1-ABzi ^ ' 

This allows us to get the relation between r and zi: 

2cosh(2r) = (1 + \A\')\z,\' + ^^^^ , (6.3) 

1 1 — /I ±3 p I Zl P 

where we have used the relation between r and the holomorphic coordinates written in (14. 2p . 
From (16. 3 p we can compute the minimum distance Vg that this embedding reaches, namely: 



cosh(2r,)^ V^VT+W , (,.„ 

Notice that this minimum distance depends on the modulus of the constants A and B , as 
well as on the phase of A B. In order to compute the function S{r) for these embeddings, let 
us compute the pullback of J. It is quite useful to work with the complex coordinates Zi and 
to obtain first the pullback of the rji forms. Actually, the puUbacks of the S0{4) invariant 
(l,l)-forms can be cast nicely as: 



I* [t]. 



2 



2 V'~ ' '""' "'"" ' \l-AB\'^\zi\y |2i|2 ' ' |2i|2 

1 ((1 + l^n - (1 + l^nil -AB\^\zi\^) dzi A dzi _ 
4 \l-AB\^\zi\^ |zi|2 

dzi A d^i 



(cosh2(2r,) - cosh2(2r)) 



kiP 



\A + B\'^ d^riAd^i ,0/ , d^i A dzi 

. PI2 — I 12 = cosh {2rg) — --^ 

ABy \zir \zir 



^* iVs) = , Vol. . r = cosh^(2r.) ^V^ • (6-5) 
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From these pullbacks we can readily compute the puUback of J, namely: 



e 2 



I* {J) 



2fc cosh(4r) + 1-2 cosh^(2r,) 29 2 cosh^(2rg) - 2\ i dzi ^ dzi 

e -TT. ^ he . - o,^ . 1^ i^^j^ 



sinh (2r) 



sinh (2r) 



Fil 



(6.6) 



Magically, the pullback of J does not contain the function a, and it is ready for comparison 
with the smeared action. In order to obtain the actual value of S{r) we need to express 
dzi A d^i = dr A d(angular). With this purpose in mind we will parameterize Zi as: 



Zi = ue 



Then, one has: 



d^i A dzi 



Fil 



and since from (16.31) it follows that: 



du 
u 



-2^ ^ A d^ 
u 



sinh(2r) 



cosh (2r) — cosh {2rq 



(6.7) 
(6.8) 

(6.9) 



we can write: 



i'{J) = 27r /dre- \e 



,^ „2fc 



IC2 J \ V2sinh(2r) 

Thus, the function S{r) in this case is given by: 



cosh(4r) — cosh(4rq) ^ a/2 tanh(2r) cosh(2r) sinh^(2rg 



sinh^(2r) ■\/cosh(4r) — cosh(4rq) 

(6.10) 



S 



^2k 



A/cosh(4r) — cosh(4rq 
V2sinh{2r) 



+ e^3 tanh(2r) 



v/2cosh(2r) sinh^(2r 



sinh (2r) A/cosh(4r) — cosh(4rq 



• (6.11) 



Plugging this result in the right-hand side of (IS.lOp and taking into account that the coeffi- 
cients of e^'^ and e^^ tanh 2r are related by a derivative: 



A/cosh(4r) — cosh(4r 



v/2cosh(2r) sinh^(2r 



sinh (2r) A/cosh(4r) — cosh(4r 



(6.12) 



one immediately gets that the profile function S{r) for this family is given by: 



iycosh(4r) — cosh(4r 



V2sinh(2r 



Qir — r„ 



sinh (2rg) 
sinh^(2r) 



(6.13) 



where we have taken into account that r > r^ on the cycle. Notice that S{r) — )• 1 as r — )■ 00 
and the massive solution becomes the solution of [10] in the far UV, as it should (see figure 
[1]). Notice also that S{r) = 1 in (16.131) for the massless case r^ = and, therefore, we recover 
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the results of [TO] in this case. As mentioned above, in appendix Owe have checked the form 
of fl and the expression of 5* for these embeddings by means of a microscopic calculation in 
the UV, where the unflavored model reduces to the "abelian" model discussed at the end of 
section [2l 

Near r = Vq the profile S{r) in (16.131) vanishes as: 



S{r) ~ 2Wcosh(2rg) ^/F^Vg , (6.14) 

which means that S{r) is continuous at r = Vg. However S'{r) diverges as 1/^r — Vg when 
r — )■ Tg. Since S'{r) enters into the energy momentum tensor of the branes (see eq. (IB. 131) 
in appendix |B]), it follows from Einstein's equations that this divergence will induce the 
divergence of the Ricci tensor at r = r^. This divergence is due to the hard- wall effect that 
we are introducing in our configuration when the flavor branes are added and it should be 
thought as the gravitational analogue of the threshold effects of field theory. In the next 
section we will propose a way to resolve this singularity in our string duals. 

7 Removing the threshold singularity 

Let us consider the class of embeddings studied in section [6l We will show how one can 
engineer a brane setup such that the unwanted singularity of 5"(r) at r = r^ disappears. 
The idea is to consider branes whose tips reach different radial positions and perform an 
average over the value Vg of the radial coordinate of the tip of the flavor branes. Actually, 
this is the way in which the threshold singularity is removed in the Klebanov-Strassler model 
with massive flavors studied in [23] . Indeed, in appendix |D] we reconsider this last model 
with the tools developed here and we explicitly show how averaging over a certain phase is 
equivalent to a particular superposition of flavor branes ending on different radial positions. 
Moreover, in appendix [D] we will also compute the function S{r) for a set of branes ending 
on a fixed Vg [i.e. the analogue of (I6.13P for the Klebanov-Strassler model) and we will 
also explore some other possibilities to perform the superposition of branes with different 
positions of their tips. 

Inspired by the resolution of the threshold singularity in the Klebanov-Strassler model, we 
will consider a flavor brane distribution containing branes with different r^'s. Furthermore, 
we will allow Vg to vary in a certain finite interval and we will weight the different values of 
rq with a non-negative measure function p{rq), which should be conveniently normalized. In 
this way the hard wall at r = r^ will be substituted by a shell of non- vanishing width. If the 
resulting profile function S and its first radial derivative are continuous the geometry will 
be free of threshold singularities. As we will see explicitly below, if the measure function is 
smooth enough the resulting profile will fulfill the conditions to have a regular supergravity 
solution. 
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For convenience let us redefine the radial coordinate as: 

X = cosh(4r) , X > 1 . (7.1) 

We will also denote Xq = cosh(4rg). We will consider distributions of branes having Xg's in 
the interval xq < Xq < xq + 6, which correspond to having quarks of different masses. The 
resulting charge density distribution is additive and can be obtained by integrating over Xq 
the profile functions f l6.13p multiplied by the measure p{xq). Since the branes with a given 
Xq only contribute to the charge density distribution S{x) for x > Xq, one has: 



Six) = I dXq p{Xq) y^^^ . (7.2) 

JXQ a/X - 1 

The measure function p{xq) must obey the normalization condition: 

/oo 
dXq p{Xq) = 1 . (7.3) 

When the measure p is a (5-function of the type p{xq) = 6{xq — Xq) the profile (17. 2p reduces to 
(I6.13P which, as we have seen, leads to background with a threshold singularity. To resolve 
this singularity we just consider measures with a finite width S and we will regard 5 as a 
regularization parameter of the threshold effect. As 5 —t- we will recover (I6.13p . Below 
we work out two simple prescriptions for the functional form of p. In both cases p{xq) is 
non-vanishing only in a finite interval xq < Xq < xq + S and the resulting S{x) and S'{x) 
are continuous and thus they source a regular geometry. Moreover, the profile functions for 
both measures are actually very similar if one compares distributions with the same width, 
as one can appreciate in the right plot of figure [H 

7.1 Flat measure 

As first example of weighting measure we consider the situation in which all the embeddings 
with different tips in the interval Xq < Xq < Xq + S weight the same. This election corresponds 
to choosing a rectangular step function in the interval Xo < Xq < xo + S which, conveniently 



normalized, reads: 



e{Xq - Xo) - Q{Xq - Xo - S) 

p[Xq) = . (7.4) 



Performing the integral (17. 2p for this measure, we get: 

2 {x- xof/ 

3 ~6V^^ 



2 jx-xpf/' ^ ^ ^A 

S[x) = t: 1^=^^ when xq < x < xq + o 



(7.5) 



2 (x - xof/' -{x-XQ- 5fl^ 

b[x) = - -^=^ when x > Xq + o 

3 5 V a; — 1 
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Figure 1: We show plots of the function S for the flat measure on the left. The red curve is 
the singular profile, the blue one is for 6 = 0.15 and the purple one is for 6 = 0.4. On the 
right, we plot, for 6 = 0.5 and 6 = 1, S for the fiat measure (red) and the peaked one (blue) 
and we see that there is almost no difference. 

and it is understood that S{x) = for x < xq. In figure [H (left) we have plotted the function 
S{x) for different values of the width 6. As shown in this figure, when 6 is increased S{x) 
grows slower in the transition region and, thus, S{x) is a milder function of x. 

Let us now consider the issue of the regularity of S{x). The potentially dangerous points 
are x = Xq, Xo + 6, where the measure p is discontinuous. It can be straightforwardly checked 
that S and its first derivative are continuous at these two points. Actually, one has: 



S{xo) = , 



S'{xo) = - 



S{xo + 5) = - 



V6 



S 



S'{xo + 5) 



3 Vx(7 

26 + 3xo - 3 
3V6{xo + 6- 1)3/2 



(7.6) 



Moreover, it follows from 07. 5p that S{x) vanishes as (x — Xq)^^^ as we approach the endpoint 
of the charge distribution at x = Xq. Thus, this profile function gives rise to a solution 
without threshold singularities, as claimed. 

7.2 Peaked measure 

In our previous example we have considered a weighting function p which is discontinuous 
at Xq and Xq + 6. We now want to explore the possibility of having a measure which 
vanishes continuously at these endpoints. To choose this new measure, we think of taking a 
distribution that reproduces a mass peak with finite width for the quark. It means that we 
choose a distribution that looks like a peak of finite width, a bit like a Gaussian function, 
but we want something simpler to be able to perform the integration. For that reason, we 
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choosqj: 



p{Xq 



The integral (17. 2p now gives: 
16 1 



Xq - a;o)(a;o + 6 - Xq) Q{xq - xq) - Q{xq -xo-6) 



(7.7) 



Six) 



157r(53/2^/^^y 



2{x^ + xl + 5^ + xqS - 2xxo - x5)E ( 



X — Xq 



\ 5 



Six) 



16 a/x — Xq 

157r52 v/^^T 



+{xo + 6- x){x -xo- 26)K ( ^^ ] 
2{x'^ + xl + 5'^ + xo6 - 2xxo - x6)E 
+ {x -Xo- S){2xo + 6 - 2x)K 



when Xq < X < Xq + 6 , 



6 



X — Xq 

6 



X — Xq 



when X > Xq + 5 

(7.8 



where, again, it is understood that S{x) = for x < Xq. The functions K and E in (17. 8p are 
complete elliptic integrals of the first and second kind respectively. They are defined as: 



E{k) 



7r/2 



dt\ 1-k sm\t) 



K{k) 



■k/2 



dt. 



1 



'0 ■ JO \/l — k sin^(t) 

If we now look at the properties of our solution at x = xq, xq + 5, we get: 

32\/5 



(7.9) 



S'(xo) = 



S"(xo) = 



S'(xo + S) 



S\xo + 6) 



IStta/xo + 5 — 1 
8(3(5 + 5xo - 5) 



(7.10) 



157rv/(xo + 5 - 1)35 ' 

which shows that S and S' are regular and, therefore, this measure leads to another solution 
free of threshold singularities. In order to compare the profile (17. 8p with the one obtained 
with the flat measure (eq. (17. 5p ) we have plotted in figured] (right) the two functions S{x) 
for two values of 5. It is rather clear from this figure that the election of p does not influence 
much S (and even less the functions of the ansatz) and, therefore, the physically relevant 
parameter is the width. For this reason, in the numerical calculations of this paper we will 
use the simpler result (17. 5p . 



8 Solutions of the master equation 

After discussing the details of our setup in the previous sections, we now move on to the 
task of finding explicit type IIB supergravity solutions. As we argued in section |3l it is 



^We thank Angel Paredes for discussions on this point. 
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enough to solve the master equation (I3.27p . since all the other functions of our ansatz will 
follow. For each P, we will have a background preserving four supersymmetries that solves 
the equations of motion of type IIB supergravity (see appendix [B]) . 

The master equation involves the profile S{r) and the function Q{r) (see (13.261) ). Notice 
that in the cases S = and S = 1, this master equation has been extensively studied in 
the literature [13]. These cases are precisely the IR and UV limits of our profiles S{r), so 
then the asymptotics of our solutions are already known. What we have to find is a smooth 
matching between them. 

We cannot provide an exact analytical solution of the master equation, but we can give 
analytical expansions in the relevant regions (around r = 0, r = tq, and r = oo), and solve 
numerically in between them. 

8.1 Analytical matching 

On general grounds we expect S{r) to be null up to a certain point r = Tq, where we have 
enough energy to start seeing the effects of virtual quarks running in the loops. Then it 
starts growing because as the energy increases it is easier to produce the quarks. It should 
eventually stabilize around S{r) = 1 since we will have enough energy so that the flavors 
appear to be massless. Although we know the specific functional form of S{r) in some cases, 
let us keep the discussion more general and assume that S{r) can be expanded in a kind of 
power series around tq as the one below: 

S(r) = 0(r - ro) [S^{r - roY/' + S^ir - Tq) + ^3(r - ro)^/^ + O {{r - rof)] . (8.1) 

It is important to notice that, according to this expansion, although S{r) is continuous, the 
[^^]-th derivative will not be, when Sn is the first non-zero coefficient of the expansion. 
Note that 5* calculated with both the fiat measure and the peaked measure are included in 
this expansion (we only have the odd coefficients for the former, and the even ones for the 
latter). 

Of course up to r = ro the solution of the master equation will be the unfiavored one. 
This solution was written close to r = in IT3l: 



P,nfl(r) = 2iVe/3 



15 V /3V 525 V 3/32 3/34 i ^^y ) 



•2) 



whercl /? > 1- Unfortunately, far from this point, it is only known numerically (except for 
the case /3 = 1, where the previous expansion truncates to the exact solution P = 2Ncr of 



^The parameter hi found in 13 is related to /3 hy hi = 2NcP. 
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[1]) . The different functions of this solution behave near r = as: 



^2h 



Nr 



/3 r^ + — (-12/32 ^ ^5^ _ 8) ^4 ^ (^^6) 



^4(#-<J.o) _ 



/3 + TT^ (6/32 -5/3-1) r2 + 



16 



15/3 



1575/3- 



■ (3/3^ + 35/33 _ 3g^2 _ 2) ^4^(^ ^^6^ 



/^ + ^(/5^-l)-^ + ^(3/3^-/^^-2)r^ + (^M 



(8.3) 



iV3/33 



1 + 



16 

9^ 



rr^ + 



32 



405/3^ 



(-15/32 + 31) r^ + O (r*^) 



a = 1 + f -2 + ^ j r^ + ^ (75/3=^ - 116/3^ + 40/3 + 8) r^ + C (r*^) . 
This solution is regular in the IR. Near r = the different curvature invariants are: 



3 Ar,'/8/32i/8 



+ 0(r) 



31p~*o 2^/2 



"'/j;/-' 



27 ivi'/^/32i/4 



^.4) 



i?,...i?-- ^ e:!:976-3072£+3456^ ^ ^^^^ 



''fiupa-' 



45 2V2Arf/^/32i/4 



From To on, 5 7^ 0, and we will have to solve the master equation with initial conditions 
given by the unflavored solution: P(ro) = Pundiro), P'(ro) = P^nfii^o)- The form of the 
solution will depend on the form of S{r) around r = Tq. 

To solve the master equation in power series close to the matching point r = tq, we need 
to know the expression for Q{r), which can be obtained from f l3.26p : 



Q(r) = N, 



2rocoth(2ro) -1 + 



sinh(4ro) — 4ro , , 2A^/- ^ ,, ,, ,n/r, 



4 (2r„ coth(2r„) - 1) N, ^^,^^pr„)S, ) (r - r^f + O ((r - r^f") 



sinh (2ro) 



2iV, 



^.5) 

To arrive at f l8.5p we have fixed the integration constant go in fl3.26p to match the unflavored 
solution at r = tq. This matching is achieved if one takes go = 0. Notice that, in Q, Nf 
appears only through the combination Nf/N^ and N^ is just an overall factor. Actually, the 
master equation fl3.27p can be written in terms of P/Nc, Q/N^, Nf/N^, and S and no other 
term depends on Nc. 
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As no term will be singular in the master equation at r = tq, the uniqueness and existence 
theorem for ordinary differential equations guarantees the existence of a unique smooth 
solution (actually as smooth as J drS) for this second order differential equation. Therefore, 
let us propose an expansion for P{r) as: 

N;'P{r) =iV-ip,„fl(ro) + N-'PU{ro){r - Tq) + Psir - rof^' + P,{r - rof+ 

(8.6) 
+ P,ir-ro)'/^ + 0{{r~ro)^) . 

Plugging the expansions (18. ip . (18. 5p and (18. 6p in the master equation (13.271) . we obtain the 
following solution: 

P^ = \ (N-'PUro) - ^ ^2) , (8.7) 

2A^/o , 8^/0. .-In/ r^^ iV-ip,„fl(ro)-2ro + tanh(2ro) 



-S, + -^S,N-'PU{ro 



^Nc 15iV, '{2roCoth{2ro)-lY-N~^Pl,{ro) 

An important lesson to extract from here is the following: our background will present 
no curvature discontinuity as long as P" is continuous (if only P' is continuous, then the 
Ricci scalar will have a finite jump at r = tq). So in this case, no curvature singularity will 
amount to having S*! = 5*2 = 0. 

In the UV we have 5* — )■ 1. The asymptotic value S = 1 will be reached exponentially, in 
a fashion that depends on the particular details of the measure used to compute S, although 
the first subleading term is universal (given by the abelian limit): 



/^+2|sh4ro _ \ ^_,. ^ ^ ^^_3.) . 



S= 1 

As we mentioned, the case S = 1 has been studied already in [13], where two possible 
UV analytical expansions were found, dubbed Class I (linearly growing P) and Class II 
(exponentially growing P). We will also have two possible UV behaviors, and the analytical 
expansions will have the same coefficients as those in [13] for the solutions with linearly 
growing P, and the same leading coefficients for the solutions with exponentially growing P. 
As argued in the next subsection, we are interested in the Class I behavior. 

8.2 Numerical matching 

If we solve the master equation (I3.27P numerically, we find, regardless of the specific profile 
S{r) we use, two qualitatively different behaviors as we go to r — )■ 00, which are in corre- 
spondence with the two classes of UV described in section 4 of [I3] . Indeed, in the deep UV, 
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Figure 2: Numerical solutions for N~^P for different values of Nf/Nc, keeping fixed the 
profile (flat measure) and tq and 6 (in the plot, cosh(4ri) = cosh(4ro) + S). The blue 
dotted line corresponds to Nj/N^ = 1. The purple line corresponds to the conformal case 
Nf/Nc = 2. And the olive dotted line corresponds to Nf/N^ = 3. Notice the expected UV 
asymptotic behaviors. 

the massive flavors we are introducing can be considered massless. We have checked that 
our numerical solutions comply with the UV asymptotic behaviors described in [13]. 

We find that in general the flavored solution only matches nicely (meaning that the solution 
will reach infinity) with the unflavored solution (18. 2 p if we choose /3 to be bigger than some 
critical value Pc, which is only known numerically (see figure [3l) and bigger than 1. This 
means in particular that the unflavored solution cannot be that of [3]. We observe the 
following: 

Assume the unflavored P up to tq is given by the numerical solution characterized in the 
IR by (18. 2p . Then there exists a (3^ such that: 

• For P < f3c, P will eventually start decreasing, crossing Q at some finite value of the 
radial coordinate and making e^^ = at that point. This solution is then singular. 

• For P = Pci P will reach infinity linearly. This solution has precisely the same asymp- 
totics as those described as Class I in [13], characterized by a linearly growing P and 
a linearly growing dilaton. 

• For /3 > Pc, P will reach infinity exponentially. This solution possesses the asymptotics 
dubbed as Class II in the previous reference, characterized by an exponentially growing 
P, and an asymptotically constant dilaton. 
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Figure 3: We plot the different values of /3c — 1 as one varies the ratio Nf/N^. The different 
curves are for different quark masses: moving from the upper curves to the lower ones, the 
values used are tq = 0,0.15,0.3,0.7,1.2; and fixed width 6 = 0.2. Notice that as the r^ 
increases (the mass increases), the growth of Pc with Nf/N^ is less and less noticeable, and 
the solution in the unflavored region is almost that of [4J (/3c ~ 1). This was to be expected 
since the more massive the flavors, the less they affect the IR dynamics. 



So the IR expansion (18.21) can be connected with any of the two known UV behaviors 
as long as we choose the parameter /3 appropriately. For an interpretation of our solutions 
as gravity duals of A/" = 1 SQCD we are interested in the ones with asymptotically linear 
dilaton [10], i.e. the ones which have /3 = /3c. Notice that the IR effects of the flavors will 
be codified in the dependence of /3c with Nf/Nc. We can then regard /3c as a measure of the 
deformation induced by the flavors in the IR. In figure [3] we explore the dependence of /3c on 
the number of flavors and their mass. 

Even if we fix /3 = /3c, and for a given ratio Nf/N^, we can still play with several parameters 
in the profile S{r), like tq, S or even with the functional form of S itself. The reader may 
wonder what would be the effect of that. We find that the qualitative behavior of the metric 
functions does not change. For instance, varying the width of the mass distribution of the 
quarks 6, just makes more or less sharp the transition from the unflavored region to the 
flavored one. We gathered in figure H] the plots of the various metric functions for some 
particular values of the parameters, just to exhibit explicitly this transition from unflavored 
to flavored background that happens around tq. 
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Figure 4: Metric functions for a case with Nf ^ 2Nc. We have used the flat measure profile 
with To = 0.5, 6 = 0.5. All the functions have the expected asymptotics. Notice in particular 
the linearly growing dilaton, in red. 



8.3 The solution for massless flavors 

Let us take rg — >■ in our expressions, keeping a finite width 6 for the measure (recall that 
also taking 5 — )■ gives back the singular solution of |lOj). This makes the lightest quark we 
are introducing massless. Nonetheless, due to the non-zero width, some of the quarks are 
massive; notice however that their mass can be chosen to be as small as one wants. In that 
respect, this solution is not a typical massless-flavor solution, as in [TOj . 
Let us consider the following expansion for the profile function S{r): 



S{r) = Sir + S2r^ + S3r^ + 0{r^ 



^.9) 



We set the first coefficient to zero because we are imposing S'(O) = 0. This expansion 
encompasses the results coming from the two measures we chose in section [71 

For this S{t\ we have to integrate the differential equation for Q in a series expansion. 
We get: 



Q[r) = N, 



A 2 NfS, 3 
3 2Nr 



16 2NfS2 
45 ^ 5Nr 



r^ + 



2Nf Si Nf Ss 
9Nr 3N7 



r^ + 0{r^ 



(8.10) 
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The expansion we get for P now is: 



N. 
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where /3 is a free parameter. We find the following IR asymptotics for the metric functions 
and the dilaton: 
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(8.12) 

Solving the master equation numerically, we find the same UV behaviors as in the previous 
subsection, that is, P grows either linearly or exponentially as r — t- 00. Again, the linear 
behavior can only be reached by choosing /3 equal to a critical value (3c (see figure E]). 

We have checked that the solution above presents no curvature singularity in the IR if we 
choose Si = 0. For instance, the Ricci scalar near r = is given by: 



R = 3e-*«/2. 



NfSi 



I 



25/4jy13/8 /513/8 r 



C(r°) 



(8.13) 



and the metric is clearly singular at r = if 5*1 7^ 0. 

Note that what is done in this subsection might be thought as a regular way to introduce 
massless flavors, as opposite to what happens in [10], where the geometry is singular in the 
far IR. This statement should be taken with a grain of salt: if we interpret Nf S{r) as giving 
the number of flavors that are effectively massless at a given scale r (see next section), we are 
clearly reading off S'(r = 0) = that there are no massless flavors in the far IR. But given 
that the tips of some branes reach the origin of the space, there are certainly massless quark 
states in the dual theory. One could conjecture about the existence of some field-theoretical 
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counterpart of the fact that the tips of the branes should be spread in order not to generate 
a curvature singularity. Unfortunately, we cannot assert any strong claim in this regard. 

9 On the dual QFT interpretation 

Up to now we have only dealt with the problem of finding a regular supersymmetric solution 
for unquenched massive quarks. This solution should capture, in a holographic setup, those 
flavor effects for which the fact that the fundamentals are massive is important. Notice 
that, in the UV, our solution reduces to the one in [10]. Thus, we expect our formalism to 
be relevant in the description of the IR physics of the model. In this section we work out 
explicitly some of the effects of massive flavors. In section 19.11 we analyze the realization of 
Seiberg duality in our massive solutions. In section 19.21 we study we study the Wilson loops 
in our background. Finally, in section 19.31 we perform the calculation of k-string tensions 
within our formalism. 

9.1 Seiberg duality 

Seiberg duality is an interesting feature of A/" = 1 four-dimensional gauge theories with 
flavors. In this section, we briefly comment on the particularities of Seiberg duality in the 
presence of massive flavors, and explain how these features are realized in our holographic 
setup. 

In his original paper [28], Seiberg argued that the IR dynamics of SQCD could be under- 
stood with the usual "electric" description, that of an SU{Nc) gauge theory with Nj flavors; 
or alternatively via a "magnetic" description, consisting of an SU{Nf — N^) gauge theory 
with Nf flavors interacting with some gauge singlets. The global anomalies of both the 
electric and magnetic theory match, a precise dictionary between gauge-invariant primary 
operators can be found (in particular the gauge singlets of the magnetic theory are related 
to the mesons of the electric theory) so that the 't Hooft anomaly matching conditions are 
satisfied, and the deformations of the two moduli spaces can be put in correspondence. 

As discussed in Seiberg's original work, to understand the effect of giving a mass to the 
fundamentals, we can just give a mass to, say, the Nf-th quark flavor. We can now think 
what happens to the IR theory both in the electric and in the magnetic picture. In the 
electric picture, the massive flavor will be integrated out in the IR, so that the effective 
electric theory will have Nc colors and Nf — 1 flavors. From the magnetic perspective, the 
mass term becomes, after working out the F-term equation of the gauge singlets, a VEV for 
the magnetic A^/-th quark. The gauge group SU{Nf — Nc) is then broken down through the 
Higgs mechanism to SU{Nf — N^ — 1), and thus the IR theory will have Nf — N^ — 1 colors 
and Nf — 1 flavors. This magnetic effective description is precisely the dual of the electric 
one just described. 
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The lesson we should extract is that Seiberg duality in the presence of massive flavors 
works very much like in the case of massless flavors, but instead of the usual duality relation 
{Nc, Nf) ^ {Nf - Nc, Nf), one should have {N^, Nf) ^ {Nf - N^, Nf), where Nf is 
the number of massless flavors. Let us see how this feature is codified in our holographic 
dual. 

The field theory whose dynamics our supergravity background is capturing is not exactly 
SQCD, but rather SQCD plus a quartic superpotential {W ~ kQQQQ). This field theory 
exhibits "exact" Seiberg duality [29], meaning that the duahty holds along the RG flo'fto and 
not just at the IR fixed point, and the theory is Seiberg self-dual. The way this was seen in 
the solution with massless flavors [10] (see [I2] for a deeper and more subtle analysis) was 
by realizing that the BPS system is invariant under the change (A^c, Nf) ^ {Nf — Nc, Nf). 
It corresponds to Q{r) ^ —Q{f) which leaves the master equation invariant. That is, the 
same supergravity solution [i.e: the same physics) can have two different interpretations, 
one electric and another magnetic. 

In the solution for massive flavors of the present paper, it can be easily seen that the 
master equation fl3.27p is invariant under {Nc, NfS{r)) ^ {NfS{r) — A^^, NfS{r)) (recall 
equation (I3.26P ). Taking into account that NfS{r) is precisely counting how many flavors 
are effectively massless at a given energy scale, this is exactly what we were expecting to 
find from the discussion above. Note that the change {Nc, Nf) ^ {Nf — Nc, Nf) is NOT a 
symmetry of the master equation with massive flavors. 

9.2 Wilson loops 

We would like to look now at the behavior of the quark-antiquark potential in the field 
theory dual to our supergravity solution, that can be studied within the gauge/gravity 
correspondence. This topic has been treated already in [10] for the case of massless flavors, 
and extended in [30] to the case of massive flavors. The interest of revisiting this calculation 
is the following: 

In the original case of [10|, the metric presented a curvature singularity in the IR; as shown 
in [31], this invalidates the gravity calculation at low energy. One could think of giving a 
mass to the quarks, that will remove the aforementioned IR singularity, but a proper gravity 
model for the theory with massive flavors was not available. The authors in [30] proposed to 
use a gravity solution built out of a flavorless solution and a solution with massless flavors, 
glued at some finite r = Vq, modeling a mass rriq ~ r^ for the quarks. This solution would 
correspond to taking in our formalism S = Q {r — r^), that introduces a very ugly curvature 
singularity at r = r^. 



■^The idea is that the QQQQ term of the electric theory becomes a mass term for the gauge singlets of the 
niagnetic theory. They can be therefore integrated out, leaving a SQCD theory with a quartic superpotential 
as the electric one. 



35 



With our gravity solution at hand, we can address the study of the quark-antiquark 
potential in a singularity- free context. Before going on, let us state that the results we 
obtain are in qualitative agreement with those of [30] , where they found that the "connected 
part" of the static potential between two non- dynamical quarks (i.e: without taking into 
account the decay into mesons) went from a Coulomb-like law at short separation distances 
to a confining behavior in the IR. Moreover, depending on the mass of the quarks niq, there 
was a first-order phase transition between these two different behaviors for masses below a 
certain critical mass rric. 

The quark-antiquark potential can be extracted from the expectation value of a Wilson 
loop, and the procedure for computing the latter within the gauge/gravity correspondence 
is well known. The idea is to introduce a probe flavor brane at r = oo (so that the probe 
quarks have infinite mass and are non-dynamical) extended along the Minkowski directions 
as well as wrapping a certain two-cycle in the internal manifold. We attach then a string 
to this brane, that will hang into the ten-dimensional geometry, reaching a minimum radial 
distance tq. We have to compute the energy E of the string and the separation L of the 
quarks at the end-points of the string for different tq's. We briefly summarize the relevant 
formulae. For details one can have a look at [31]. Defining: 



where C = f{ro) and we are using string frame, we have that 



(9.1) 



Attaching a string to the probe flavor brane we are introducing can be done whenever it 
is possible to impose Dirichlet conditions on the string end-points. For our geometry, as 
discussed in [31], this is possible when lim V{r) = oo. Since for large r, V^ ~ e*"'^, this 
conditions holds only for the solutions with an asymptotic linear dilaton. For these solutions 
we plot the results in figure O 

As mentioned above, the quark-antiquark potential exhibits two different behaviors: an 
inverse-power law in the UV, and confining in the IR, where the massive quarks have been 
integrated out, and the dynamics of the unflavored theory is recovered. The transition 
between these two behaviors can be smooth as in the plot on the right, or a first-order phase 
transition (the derivative of the energy has a finite jump), as in the plot on the left. As 
explained in [30|, this behavior could be expected whenever we have two scales in the theory. 
In the present case, these scales are the gaugino condensate and the mass of the quarks. 
More precisely, our background does not have a sharp value for the mass of the quarks, but 
rather a distribution of masses with a certain width. The phase transition shows up when 
the mass of the heaviest quark is smaller than a certain critical mass rric, set by the gaugino 
condensate. 
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Figure 5: We plot the energy of the Wilson loop E vs. the quark separation L. We have fixed 
Nf/Nc = 1 and tq = 0.05. For small widths of the brane distribution we have the plot on 
the left, where we observe a first order phase transition. As we increase the width, the mass 
of the heaviest quark becomes of the order of Asqcd, and the phase transition disappears, 
as shown in the plot on the right. Notice the similarity between these curves and the G-P 
(Gibbs free energy vs. pressure) curves of the Van der Waals gas. 



We could have pursued a more detailed study of these phenomena, like analyzing the de- 
pendence of rric with the number of fiavors, their masses, and their distribution, or exploring 
the decay into mesons characterized by the string breaking length. A fast analysis revealed 
that the way we distribute the quarks is not very relevant for these observables, and that 
other features follow qualitatively the behavior described in |30]. We would like to stress 
though, that the present calculations are performed in a background without any pathology, 
giving a more solid foundation to them. 



9.3 k-string tensions 

One of the most interesting features of the IR physics of the confining A/" = 1 theories is 
the existence of the so-called k-string states, i.e. of fiux tubes induced by sources with k 
fundamental indices. It was argued in [32] that such a state can be described by D3-branes 
extended along one of the Minkowski spatial directions, time and wrapping a two-sphere 
in the IR geometry. For the unfiavored geometry of section [21 the tensions of these k- 
strings obey a sine law. It is important to notice that, in order to get the results of [32] . 
it is crucial to find the RR two-form potential C(2). In the approach of [32] the potential 
C(2) is converted into a NSNS two-form 5(2) by means of an S-duality transformation and, 
then, the fiux stabilization mechanism of [33] is applied to determine the configurations that 
minimize the energy and to obtain the corresponding tensions. In our fiavored background 
the Bianchi identity of F(3) is violated due to the presence of D5-brane sources (eq. (13. 7p ) 
and, accordingly, one cannot define the RR potential in regions in which Q is different from 
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zero. However, in our massive flavored case, the probe D3-brane will only explore the deep 
IR region near r = 0, where there are no flavor brane sources since the proflle function S{r) 
vanishes there. For this reason we will be able to deflne the potential C(2) in this region 
and we will proceed with the analysis of the k-string states. Notice also that, in our low- 
energy analysis, one would not expect to flnd k-string breaking due to quark-antiquark pair 
production. However, we will clearly flnd screening effects due to quark loops which will 
modify the tensions. 

Let us begin our analysis by studying the IR geometry near r = 0. Following [lOj, we will 
consider the submanifold |34] deflned by the conditions 6 = 6, (p = 2tt — (p a.t r = 0. From 
the IR behavior (18. 3p of our solutions it is straightforward to verify that the metric (13.81) 
along this submanifold, in the string frame, takes the form: 



ds^ = e*(°) 



dx?_3 + N.pUx^ + sm^x{d6^ + sm^6d(f)^) 



(9.3) 



where the angle x is related to the coordinate ip by means of the relation: x = (z/' — 7r)/2. 
Clearly, the angles x? 6 and parameterize a non-collapsing three-sphere at the origin r = 
and we should take x to vary in the range < x ^ ^- Notice also that the constant 
P characterizes the size of this three-sphere. At r = the charge density of the flavor 
branes vanishes and, as a consequence, there is no violation of the Bianchi identity of -F(3) . 
Therefore, it will be possible to represent at this point -F(3) in terms of a two-form potential 
C(2) (-^(3) = dC(2)). Actually, it is straightforward to check that C(2) at r = in these 
coordinates takes the form: 

C(2) = -N,C{x) sm6d6 A d0 , (9.4) 

with C{x) being the function: 

C(x)=-X+^. (9.5) 

Contrary to the approach followed in [32], we will perform our analysis directly in the D5- 
brane background, without performing the S-duality transformation (see also [35]). Accord- 
ingly, let us now consider a probe D3-brane moving in our background. Its dynamics would 
be governed by the action: 



Sd3 = -n /" d^e e"* ^- det {g + F) + n J F A C(2) 



(9.6) 



with g being the induced metric on the worldvolume of the D3-brane and F the worldvolume 
gauge fleld. We now consider that the D3-brane is extended in (t, x, 6, 0) in the metric (19.31) 
at r = and at flxed values of x cind of the other two Minkowski coordinates. We will 
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also assume that there exists an electric worldvolume gauge field Fq^ along the Minkowski 
direction. In this case, the D3-brane action can be written as: 

5^3 = f dtdx C , (9.7) 

where we have integrated over the angles {6, 0) and C is the effective lagrangian density, 
given by: 

C = -AtxT; N, [/3 sin2 x ^e2*(o) - F^ + Fq. C{x) ] ■ (9.8) 

The equation of motion for the electric worldvolume field is: 

dC 



dF< 



constant , (9.9) 



Ox 



which is nothing but Gauss' law. Following [36j the constant on the right-hand side of (19. 9p is 
fixed by imposing the quantization condition corresponding to having k fundamental strings 
along the x direction: 

kTf , keZ, (9.10) 



dFo, 

where Ty = l/(27ra') is the tension of the fundamental string. This condition determines 
the electric field in terms of the angle x- Indeed, let us define a new function C{x) as: 

C{X) ^ C{x) + ^ . (9.11) 

Then, one has: 

e*W C{x) 



Ox 



, , = . (9.12) 

v//32sin^X + C(x)2 

Notice that Fq^ is the momentum of a cyclic coordinate that can be eliminated from the 
lagrangian. The correct way to do this is by performing the Legendre transformation and 
computing the hamiltonian as: 

dC 



H = dx 



Ox 



- C 



(9.13) 



dFo. 

By calculating explicitly the right-hand side of this equation and writing the result in terms 
of X, we get: 

H ^ ,.ne^'^> N. I ,. ^F^;^^7^, . ,9.14) 

Let us minimize the energy with respect to x- For this purpose it is interesting to notice that 
the function C{x) satisfies dC/dx = — 2sin^x. Using this property of C it is straightforward 
to prove that, for a given integer k, the energy is minimized for the Xk which satisfies: 

C{xk) = y sin(2xfc) , (9.15) 
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or equivalently: 

Xfc - ^ + ^^ sin(2xfc) = , (9.16) 

which is the equation written in |32] with /3 instead of the b of [32] . It is also immediate to 
find the tension of the fc-string object, namely: 

e*(o) A^ , 

Tk = -ir^ /5 ^'^Xk Vl + (/32 - 1) cos2 xk . (9.17) 

The worldvolume electric field corresponding to this solution is: 



Ox 



Vl + {13' - 1) cos2 Xk ■ 



(9.18) 



It is interesting to point out that (I9.16P does not change under the transformation k — )■ N^ — k 
and Xk — ^ T^ — Xk- One can also check that the tension in (I9.17P does not change under 
this transformation, while the electric field (I9.18P changes its sign. Notice also that in the 
unflavored case reviewed in section [2] one has (3 = 1 and we recover the results in [32 | [36 | [35] . 
The case with (3^1 for the generalized unflavored models with the IR behavior (18. 3p was 
considered in ^U\. Notice that, in our case, the parameter /3 is related to the mass of the 
quarks and to the number of flavors by means of the matching conditions discussed in section 

El 

Let us look at the tension of the k-string as a function of k/N^, for different values of p. 
First, we need to solve (I9.16p . Depending on the value of (3 and k/Nc, we find that there can 
be up to three different solutions. We then have to check which one corresponds to the true 
minimum of the energy (I9.14p . We notice that for k/N^ < 1/2, the minimum of the energy 
is given by the solution for Xk closest to as we can see in figure [6l while for k/N^ > 1/2, it 
is the solution closest to vr. 

Knowing the correct value of Xk for each /3 and k/Nc, we can go on and plot the tension 
(I9.17P as a function of k/Nc for various values of (3, as shown in figure [3 When /3 is close to 
1 (which is the smallest value it can reach, corresponding to the unflavored case), the tension 
of the k-string can be approximated by: 

T, . £!i^,,.i„ (It) . ,9.19) 

27rV ^ \NcJ ^ ' 

Thus, in this low /3 case, the screening effect due to the flavor is manifested in the tensions by 
just multiplying the sine formula by the deformation parameter /3. In turn, /3 can be related 
to the number of flavors and their masses by means of the matching condition studied in 
section [8l Notice that, for a given number of strings fc, the tension of the flavored k-string 
is higher than the one corresponding to the unflavored theory. Actually, this is what is 
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Figure 6: The (lower) red curves correspond to the function on the left-hand side of f l9.16p . 
whose zeros are the angles that extremize the energy (19.141) . The (upper) blue curves are 
plots of the energy. We are taking in all cases k/N,, = 0.4. On the left, we plot these curves 
for /3 = 1.2 and we see only one zero, which corresponds to the minimum of the energy. 
On the right, for /3 = 2, we notice the appearance of two new zeros, one corresponding 
to a maximum of the energy, and the other one to a metastable configuration. The true 
minimum, however, moves towards x = 0? that is towards the north pole of the sphere, as /3 
is increased. 



expected on general grounds since the screening reduces the (negative) binding energy and, 
therefore, it increases the total energy {i.e. the tension). 

As P goes to infinity, the binding energy becomes smaller and smaller and the tension of 
a k-string is a linear function of k. In this case one can analytically obtain the approximate 
solution of (I9.16P which corresponds to the minimum of the energy. Indeed, if /3 is large 
the only possibility to solve (I9.16P is by having sm{2xk) small. One can show that when 
k/Nc < 1/2 this equation is solved for Xk ~ 7ik/(3^Nc, while for k/Nc > 1/2 the energy is 
minimized for Xfc ~ tt — 7r{Nc — k)/P'^Nc. The corresponding tensions in these two cases are: 



-m T, 



kTf 



{N,-k)Tj , 



for 0<k< NJ2 



for NJ2 <k<N^ 



(9.20) 



which shows that, when /3 is large, the screening effects are so large that the binding energy 
is very small and one can regard the flux tube as composed by non-interacting strings with 
vanishing binding energy. One can visualize this behavior as /3 is increased in figure [71 

10 Conclusions 

Let us summarize our main results. We have considered the addition of unquenched massive 
fundamental matter to the gravity dual of the TV = 1 SQCD-like theory obtained when 
D5-branes wrap a two-cycle inside a Calabi-Yau threefold. The matter fields are added by 
means of D5-branes that are extended along a non-compact two-dimensional submanifold of 
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Figure 7: This plot corresponds to the tension of the k-string as a function of k/Nc- The 
values for /3 are 1 for red, 1.1 for purple, 1.4 for blue and 2 for orange. 



the internal space. These flavor branes do not reach, in general, the origin of the holographic 
coordinate r and their charge density depends on r. In order to incorporate consistently in 
the backreacted background the effects of this dependence we have modifled in a non-trivial 
way the ansatz of [lOj, by including new terms in the RR three-form F(3) which depend 
on a proflle function S{r) and on its derivative. We have shown that the BPS system can 
be reduced to a master equation, containing S and S", which is a generalization of the one 
found in [13]. This equation can be integrated numerically and, by matching the unflavored 
solution at the scale at which S" = 0, one can find a supergravity solution in the whole range 
of the radial coordinate. 

Our solutions satisfy Einstein equations with sources and one of the non-trivial points 
we have addressed is the determination of the distributions of branes whose charge density 
and backreaction have precisely the form that we have adopted in our ansatz. We have 
verified this fact by means of a macroscopic calculation (comparing the action of the full set 
of branes with the one corresponding to a representative), as well as by a direct microscopic 
calculation of the charge density in the UV. We have also shown how to get rid of the 
curvature singularities which appear at the position of the tip of the branes. 

After all these developments we have been able to find regular supergravity backgrounds 
dual to A/" = 1 SQCD-like theories with massive unquenched fiavors. Our results generalize 
those in [10] in the sense that our solutions incorporate the effects of the mass scale introduced 
by the quark mass and, at the same time, they resolve the IR curvature singularity that limits 
the applicability of the geometry of [TO] to explore holographically the M = 1 SU{Nc) gauge 
theory with flavors. 

We have studied some observables of the fleld theory for which the IR structure and the 
mass scale introduced by the quarks are relevant (see [37] for further analysis along these 
lines). We think that our formalism provides a framework to explore holographically the 
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A/" = 1 SQCD-like theories in a firmer basis. 

For the standard field theory interpretation, we focused mainly on solutions with linear 
dilaton at infinity. However, it is also possible to get backgrounds where the dilaton is 
bounded. On such solutions, we can apply a transformation found in [SS], called rotation. 
It would give us a completely regular version of the background of [38] , and for that reason 
it could allow for a better understanding of the dual field theory. 

We could also think of introducing new scales in our supergravity solutions, by choosing 
a profile function for the fiavors created by several separated shells of branes. The function 
S{r) would then look like a series of steps, whose lengths would be related to the different 
positions of the flavor branes. If, by engineering the proflle in a smart way, we can have the 
/9-function of the theory almost zero for one of the steps, the this could lead to a dual fleld 
theory exhibiting a walking behavior [391 EB HOI SI] • 
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A Super symmetry analysis 

In this appendix we will study the realization of supersymmetry for a background of type 
IIB supergravity with metric and RR three-form F(3) as given by the ansatz written in eqs. 
f l3.8p -( l3TT3l) . To perform this analysis we will choose the following basis of vielbein one-forms 
for the metric (13. 8p : 



^^ 



/dx\ 2 = 0,1,2,3, 

e't> = e^+^ sin ^ d0 , e^ = ^(w^ + a(r)d^) , ^^'^^ 

g/+9 _ , g/+fc 

e^ = (w^ - a(r) sin 6* d0) , e^ = {co^ + cos9d(f)) . 
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In order to find the supercharges preserved by a given background of type IIB supergravity 
one should find the KiUing spinors e which leave invariant the dilatino and gravitino field 
under a supersymmetry transformation. These spinors are characterized by a set of com- 
muting projections. In the case of the unfiavored background reviewed in section [2] these 
projections were found in [8]. We will impose that the spinors e of the fiavored background 
satisfy the same type of conditions as in [8J , which can be written as: 

ie* = e, Te<f>e = Ti2e, T^me = (cosa + sinar02)e , (A. 2) 

where Ta^a2--- denotes the antisymmetrized product of constant Dirac matrices and the Oj's 
are fiat indices corresponding to the vielbein (lA.ip . In (IA.2p a = a (r) is an angle to be 
determined. The result of the analysis of the supersymmetry variations can be nicely recast 
in terms of the two fundamental forms of the underlying geometric S'[/(3)-structure of the 
internal complex manifold. These forms are the (1, 1) two-form J and the holomorphic (3, 0) 
three-form Vthou which are defined in terms of fermionic bilinears as: 

with e being a Killing spinor normalized as e^ e = 1 and e"^^"^" = e"^ A e"^ ■ ■ ■ . Actually, 
by using the projections l\A.2\} satisfied by e one can express the S'f/(3)-structure forms as: 

J = e''^ + (cos ae''' + sin a e^) A e^ + (- sin ae'^ + cos a e^) A e^ , 

^hoi = e^-^+*/^ {e'' + ie^) A ((cos ae^ + sin ae'^)+i e^) A ((- sin ae'^ + cos ae'^)+i e^) . 

(A.4) 
The conditions imposed by the preservation of A/" = 1 supersymmetry can be written as the 
following set of equations to be satisfied by the structure forms: 

dfihoi = , , , 

d(e^OAJ)=0, 

where *6 denotes the Hodge dual with respect to the internal part of the metric (13. 8p . 
Plugging in (JA.51) the explicit expressions of J and Qhoi written in eqs. (lA.4p . as well as our 
ansatz for the metric and RR three-form, one gets the following system of first-oder BPS 
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equations: 



/' 



h' 



$' 



-g-h 



4 8 sin a 



-2g-h 



4 sin a 



2Nc{b-a) + 2a e^^ cos^a + 2e^+^ sin(2a) 

((b - a)N, + (e^s + e^^)a ] (1 + cos^ a) 



I -2q 

g = e ^ 



2 e'' (2 e^^ + e^'^ - A^^) sin(2a) + 2e^ 

,9-h 



sma 



fc' 



-g-h 



sma 
2e-2^ 



(e^'^ - Nc) cos a + — - [ (a - 6)Ar, + (e^^ - e^'')a] si 

{bN, + {(?^ - N,) a) cosa + 2 e'"'^ (e^^ + e^'' - N,) sin 
2e^s 



a 



b' = —Li H — — — [a cos a + 2e^ ^ sina] , 
together with the following two algebraic relations: 



(A.6) 



cot a 



4a 



(e^^a^ -46^'^ -e^^) 



^g-Zi 



2 [2 



cos a — e5-'*asina 



'] 



oa+h 



(L2 + 2) sin a + 4—— (a cos a + 2e^ ^ sina) 



(A.7) 

By computing the derivative of the second of the constraints in (1A.7I) . one gets the remarkably 
simple relation between the derivative of L2 and Li. 



L'2 = Ue^-3 coia - 2a) Li 



1 + a2 + 4 e2^-29 



Li, 



{A.t 



where in the last step we have used the first equation in flA.7l) to eliminate a. 

A.l Partial integration 

Although we have not been able to get an analytic expression of the general solution of the 
system (1A.6|) . we have been able to simplify it and to perform a partial integration. Indeed, 
let us proceed as in [10] and begin by defining two new functions C{r) and S{r) as: 



C 



s 



1 + a2 + Ae'^'^-^a 
2a 



+ 2a2(4e2'^-29-l) + [i + e^f'-^ay 
2a 



(A.9) 
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From the BPS equations (1A.6P and (\A.7\i one can show that C and S satisfy the following 
second-order differential equations: 



C" - 4C = , S" - 4S = , C^ - S^ 

These equations can be immediately integrated as: 

C = cosh(2r) , S = sinh(2r) , 



1 . 



(A.IO) 



(A.ll) 



where we have fixed the integration constants by imposing the r — )■ behaviors C* ~ 1 + 2r^ 
5" ~ 2r. Using the definitions of C and 5* in flA.9|) . we get the following relations: 

a^ + 1 



. ^29 ^ g2h 



- e^^ cosh(2r) 



a' + 2a' {Ae'^-'^ - l) + (l + e''"'^)' = ^a' sinh2(2r) . 
The relations (JA.121) can be combined with the first constraint in ( ]A.7l) to give: 

a'-l 



=29 



e'^ = e 



^^3 a cot a . 



Indeed, by summing this last equation and the first one in (lA.12p . we get: 

p9-h 



cot a 



(a — cosh(2r)) . 



(A.12) 



(A.13) 



(A.14) 



Using this result it is easy to find the following expressions of sin a and cos a in terms of the 

other functions: 

2e^-3 cosh(2r) -a , , . 

sin a = — — - , cos a = — — r — . (A. 15) 

sinh(2r) sinh(2r) 

The following useful combinations of cos a and sin a can be found from the previous equa- 
tions: 

1 — a cosh(2r) 



a cos a 



+ 2e^-^ si 



sma 



sinh(2r) 

2 cos a — e^~ a sin a = 2 coth(2r) . 
Using these results, the equation giving h' in the system flA.6P reduces to: 

, _ 2e'3 1 _ a cosh(2r) 
~ ^Vc sinh(2r) ^ ' 

while the second equation in (IA.7P becomes: 



(A.16) 



(A.17) 



1 



2cosh(2r) 



26^3 

L2 + 2 + -— (acosh(2r)-l) 



(A.18) 
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Moreover, the relation (lA.Sp takes the form: 

L2 = -2cosh(2r)Li 



(A.19) 



which is nothing but fl3.15p . By combining ( ]A.17p and (lA.lSP we can obtain the following 
differential equation for b: 



, , , 2 cosh(2r) , L2 + 2 

+ = 

sinh(2r) sinh(2r) 

Using ( jA.igp this equation reduces to: 



L2 + 2 



2cosh(2r), 

b -\ ^ — -b 

sinh(2r) sinh(2r) 



+ 



Li 



a 



2cosh(2r) 



(A.20) 



(A.21) 



Applying the method of variation of constants, we can integrate (IA.21I) and get 6 as a function 

of L2. The result is: 

2r + T]{r) 



smh(2rj 
where T]{r) is defined as the following integral involving L2: 

tanh(2p) 



ri{r) 



dp 



L2{p) + 



L^ip) 



Performing a partial integration, one can rewrite ri{r) as follows: 

1 



7][r) 



tanh(2r)L2(r) + / dp tanh^(2p) L2(p) 



(A.22) 



(A.23) 



(A.24) 



Plugging in flA.24p the definition fl3.16p of the profile function S we just get f l3.20p . We can 
check that the above formulas give the right results in the known cases. For example, in the 
unflavored case L2 is zero and it is obvious from (1A.24P that rj vanishes and one recovers the 
correct result of section |2l Moreover, in the massless case studied in [10] , the function L2 is 
constant. Actually, in this last case L2 = —x with x = Nf/N^.- From flA.23P or flA.24p we 
get rj = —xr and the result of eq. (4.14) in [10] is recovered. 

There are combinations of the BPS equations in (1A.6P that become particularly simple. 
One of these combinations is: 

2$' + h' + g' + k' = 2coth(2r) . (A.25) 

Integrating this relation, one gets eq. f l3.22p . Other interesting combinations are: 

^^'^ - Nc) cosh(2r) - ae^^ + N^b 



dr sinh(2r) 



d 

dr 



(A.26) 



[ae'^] 



sinh(2r) 



e2fc - A^^ + e^^ (1 - acosh(2r)) + N^b cosh(2r 
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Combining the second of these equations with (lA.lSp . one can immediately compute the 
derivative of the function P defined in (13.231) . with the result: 



P' = 2 e'" - NfS 



(A.27) 



which is just the second equation in fl3.25p . Moreover one can also compute from flA.261) and 
( JA.181) the derivative of Q. The result can be written as: 



d 



Q 



dr Lcoth(2r) J coth^(2 



2N^ -NfS 



(A.28) 



whose integration yields fl3.26p . Another combination of derivatives that becomes particu- 
larly simple when one uses the BPS equations is the following: 



h' + g' — k' = ae 
From this equation we easily get: 



2k-2h 



2coth(2r) . 



1 



dr 



4sinh^(2r 
which is nothing but the master equation fl3.27p . 



sinh (2r) — 



P'-Q' 



P' + NfS 



(A.29) 



(A.30) 



B Equations of motion 

In this appendix, we state for completeness the equations of motion. Let us first rewrite the 
action, constituted of the type IIB action and the source action: 



•S* — Sub + S. 



IIB ~r »Jsources 



where: 



and: 



S 



IIB 



2^10 



Iv^{r- ld,<!>d^<!>^ -\j (e*F(3) A *F(3)) 



(B.l) 



(B.2) 



^sources = -T5 j {c^^lC - ^(6)) A fi , (B.3) 

with K, being the calibration form for the D5-branes which is related to the 5'f/(3)-structure 
as: 

/C = e*d^x A J . (B.4) 

First we give the modified Bianchi identities for the flux: 



dF(3) = 47r2fi 



(B.5) 
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where we used that 2k1qT^ = An"^. The equation of motion for the flux reads: 

d (e* * F(3)) = . (B.6) 

For the dilaton and the Einstein equations, we define first the following notation: 

^{py\p) = ^^'^"••'^^A^,...^^ , (B.7) 

for any two p-fornis W(p) and A(p). One can easily prove that, in a ten-dimensional manifold, 
one has: 

/ W(p) A A(io-p) = - / ^/^\J{*UJ) . (B.8) 

Using these results, we can write the equation of motion of the dilaton as: 

-^d, {./^g>^-^d,<^) = i-e*F(|) - 27r2e*/2fiX*/C) . (B.9) 

Finally, the Einstein equation is: 

R,u - \g,uR = \d,^d,^ - \g,.d,^d,^ + ^e* [QF^p.F/'^ - 9,.F^,)) + T^r"'^ , (B.IO) 
where 7^°"''^'^'^ is the energy-momentum tensor coming from the source action. It is given by: 

T^T'"'' = y e*/2 (6^,.fi.(*/C) - n,,,,,,, (*/C)/^''^^^ ) . (B.ll) 

From (IB.lOp and (IB. lip , one can get an equation for the Ricci tensor as: 

- -^e*/2(2fi,p,p,p3 (*/C)/^^^''^ - 3g,M4*JC)) . 

To show more clearly the dependence of J'^°"i'ces q^ ^]^g smearing profile of the flavor branes, 
we can write T^'^^'^'^'^^ in flat components in the vielbein basis (]A.1|) : 

^sources _ ^/ -2g-2h-2k-^/2 (2k cj , 4e'^ + 6^^ (a - COsh(2r))' A 

ry-isources f_ —2g—2h—^/2 q Tisources 

-'-rr — r)^ "^"-^aS ' 

rT-isources / —2g~2k—<^ 12 qI T-isources /t3 1 qN 

smh(4r) ^^ 

THSOurces _ ^ f ^~2g~2h-2k->i>/2 ^^ + ^ (a — COSh(2rjj ^^ _ resources 

11 ~ 2^ sinh(4r) ^ " ^^ , 

resources _ ^^ f ^-g-h-2k-^/2 ^ ~ COSIl^irj ^^ ^ resources 

2 smh(4r) ^ 
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It is clear from (IB.12P and (1B.13P that to have a geometry with regular Ricci tensor one 
must require that, as stated in the main text, both S and S' are continuous functions of the 
radial variable. Moreover, it is also possible to verify that any solution of the first-order BPS 
system written in ( 1A.6I) and ( JA.71) is also a solution of the second-order equations of motion 



for supergravity plus sources written above. 



C Microscopic computation of Q 

Let us show here how the results obtained for the smearing form in section Ej using just 
the knowledge of one embedding of the family plus an ansatz for the functional form of Q 
(eq. fl3.18p ). could be derived from a purely microscopic computation: i.e., by summing the 
contributions to the smearing form of all the embeddings in a given family. 

Notice that this microscopic approach does not assume any specific ansatz for the smearing 
form. Obviously, one can expect it to be much harder to be carried out. Indeed, except for 
some very simple cases ( [2T1 [221 [23t [20] ) , a full reconstruction of the functional form of Q for 
massive quarks from the microscopic family of embeddings giving rise to it is not available 
in the literature. The use of the holomorphic structure of our internal manifold developed 
in section [H will be instrumental to carry out this microscopic computation. 

C.l Holomorphic structure in the abehan hmit 

For simplicity, we focus in this appendix on the UV limit (r — )■ oo) of our backgrounds. This 
limit corresponds to the so-called abelian solution (see the last paragraph of section [2]). The 
holomorphic structure simplifies a little bit in this limit, and one can define a new set of four 
complex variables Q (i = 1, . . . , 4) that parameterize now a singular conifold: 

ClC2-C3C4 = 0. (C.l) 

The radial variable r is related to the Q in this case as: 

Eic^i' = ^"^- (C-2) 

1=1 

The expression of these complex variables in terms of the coordinates of the internal manifold 
can be read from (14. 8p . One just needs to take the r — )■ oo limit there, to obtain: 

. r . d , 9 , i>-<t>-4> > r ^ ^ j i'+'t>+4> 

Ci = — e sm - sm - e 2 C2 = e cos - cos - e 2 

Si 2 2 ' S2 2 2 

(C.3) 

. ... 0.6 ■ 1/.+0-0 r ■ 6 ; ii-<t>+i 

G = e cos - sm - e 2 q = _e sm - cos - e 2 

^■^ 2 2 ' S4 2 2 
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This abelian geometry inherits the SU{2)l x SU{2)ji symmetry of the non-abehan one 
(actually the isometry group is enlarged to SU (2) lX SU (2) rxU (1)) . Again, taking carefulljQ 
the limit r — > oo in the non-abelian expressions fl4.13p and f l3.18p for the fundamental two- 
form J and the smearing form Q we get: 

e~~ J = —dr A (wg + cos ed(j))- — sin OdO Ad(j)- e^^ sin OdO Ad(j), 

— ^f]= sin^d^Ad^A ('Ssin^d^Ad^-S'drA Td^ + cos^d^)) . (C.4) 

Then one can define S'0(4)-invariant (l,l)-forms ?7j (i = 1, . . . ,4) as in (14. lip , and express 
both J and Q in this abelian setup as: 



e-^ J = 2ie-^'' 



e2g 



e"^ (r/i + 2 e-'^V2 - 2 e'^'^r/a) + ^ (r/i + 2 e-'^-^, + 2 e-'^vs) - e''e-'^V2 



-— fi = -8e-^'-Sr], A (r/i + ie-'^-f^^) + Se-'^^S' r,^ A (r^i - 2e-2^'r73) , (C.5) 

where the 77's are the abelian (1, 1) two-forms, which can be obtained from (I4.12p by keeping 
the leading term when r — )■ 00. 

C.2 Abelian limit of the simple class of embeddings 

Let us now calculate S{r) for the abelian version of the class of embeddings discussed in 
section [6l The first thing to notice is that the parameterization ( 16. ip is not good in the UV 
limit. Indeed, as z^ = ziZ2/zs when r — )■ 00, the two equations in (16. ip become the same. For 
this reason, to study the cycle in the UV, it is better to use instead the first two equations 
in (16. 2p and write the equation of the embedding as zi = Cz^ and Z1Z2 = fi, with C and fl 
being arbitrary complex constants. By taking the UV limit in which Zi ^ d, one concludes 
that the abelian limit of the particular representative of the embedding studied in section [6] 
is: 

(i = cc,, CiC2 = /i. (C.6) 

One nice thing of the abelian limit is that ( JC.6P can be easily solved in terms of coordinates: 



9 = 9o, = 00, and -smO e^^'e'^ = fi = -e^^'e'^ , (C.7) 

where we have parameterized the constants above as C = tan ^ e"*"^", jl = ^ (sin ^ cos y) , 
and Tg is the minimum radial distance this embedding reaches (e^'"'' = |2/i|). If we now 



*In the abelian limit: a — )■ 0, cosli2r — )■ sinh2r — >■ ^, and ae^'' — >■ 1 + 4e^'' ^s. 
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rotate this embedding with the SU{2)l x SU{2)r isometry group (see (14. 6p ). we obtain the 
expression of a generic embedding of the family as /i = and /2 = 0, with: 

_ h + aC 

Ji — ii — - — tt; Cs ; 
a — oC 

/2 = (^^\k\^-\l\^)(^(^-klC-'(^ + klCCl)- {\a\^-\b\^-abC + abC)~^fL. (C.8) 

The smearing form should be computed as an appropriately weighted sum of the transverse 
volume forms of each embedding. The formula for real constraints was first written down in 
[21] . and the generalization to complex constraints like the ones we have now is immediatqj: 

^ = T^ I dp<5^'^ (/i) 5^^^ if 2) d/i A d/i A d/2 A d/2 , (C.9) 



-2i 



c* 



where p is the (normalized to the unity) measure of SU{2)l x SU{2)b., multiplied by Nf, 
and is given by: 

dp = dadad6d6dA;dA;d/d/(5(|ap + |6p-l)5(|A;|2 + |/|2-l) ^^ , (CIO) 

A shortcut for computing flC.9P is to notice that all the embeddings of the present family, 
in virtue of the first equation in fIC.Sp . sit at constant values of 9 and (f). Since it turns 
out that the action of SU{2)l corresponds precisely to varying these constant values over a 
two-sphere, the smearing form Vt necessarily exhibits a ^ sin ^ d^ A d0 factor. We are not 
interested in getting this trivial part from (]C.9p . so we factor it out by defining an effective 



(complex) two-dimensional problem. We can define a new pair of effective complex variables: 

u ■ i>+(i) u • ^ — 4> 

^i = e"cos-e* 2 , ^2 = e"sin-e'2 , (C.ll) 

and the family of embeddings over which we want to smear recasts as 

/ = (|Ap - |i?n i^ i^ + ABe^ -ABil~p = ^, \A\^ + |5|2 = 1 . (C.12) 



(Recall that /i = ^e^^'^e*''', see (1C.7|) ). Forgetting for the moment about the correct normal- 



ization factors, the integral we want to compute is: 

W= I dAdv4d5dfi5(|A|2 + |5|2-l)5(2)(/)d/Ad/. (C.13) 

^The complex Dirac delta should be understood as (5^^''(/) ~ (5(Re(/)) (5(Ini(/)). The ^^ prefactor is 
mcluded because d/ A d/ = -2idRe(/) A dlm(/). 
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Performing this integral requires a little bit of care with the delta functions, but other than 
that, it can be considered straightforward. Let us sketch how one could proceed. To simplify 
the calculation, we reparameterize the integration variables as follows: 



A=^^^e-, 5 = ^^^e-. (C.14) 

Clearly, one has |Ap + \B\'^ = Ui and \A\'^ — \B\'^ = U2 and: 

1 /'oo /'oo /'27r /'27r 

dAdAdBdB6{\A\'^ + \B\'^-l) = - du2 dm den / das^l^i - 1) • (C.15) 

The integral in ui is then immediate. Rewriting e*"^ = X2 + i|/2, and using that: 

/ da2 = 2 dx2di/2 (^la^s + vl - 1) , (C.16) 

Jo JR2 

we can write: 

r.27r /.I 



1^=/ dai / dM2 / dx2dy2 6{xi + y^-l)5{R)6{I)df Adf, (C.17) 

where i? = Re (/|„^=i), / = Im (/|„^=i). In the new variables one has: 

-L ,u^. ^,„^ L. ^ ^o ^ ^ i .Y 



/U,=i = --e^(— ^)^1 - ul e? + «2 6 6 + -e^(-"^+"^)y^l - «1 ^2' - /^ • (CIS) 

Solving R = I = for 2:2 and t/2, the integral of the corresponding deltas produces a factor: 



dR dl 



da;2 dy; 



-1 



1 



T _ 7,2 lit 14 It |4| ' 



(C.19) 



and leaves the argument of the remaining delta function as: 

:i-«D(i6P-i6P)^ 



6{xl + yl-l) = 6i —^1±M^^^(^U2 - U2+){U2 - U2-) ) , (C.20) 



where U2± are given by: 



4|/x66|cos(^-7)±|l6P-|e2lV(|6P + l6IT-4H^ 
«2± = Tvrrr-Tr^^o • (C.2IJ 



(I6l^ + I6i 



We have at this point: 



W^ = / dai / dM2^r^PV^'^ 1 7^^i3T^^^^r^T;^^(«2 " W2+)K - «2-) ) d/Ad/ . 



(C.22) 
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In this expression, nothing depends on cti, so one can integrate it easily. Also, both M2+ 
and U2- are between —1 and 1, so they both contribute to the integral. Using (IC.2I) . and 
replacing U2+ and U2- by their values flC.21|) . we finally get: 



,,, , l-cos^ + 2e^''«-^''cos^ ,, ,- , l + cos^-2e^''«-^^cos^ ,, ,- 
W = An , d^i A d^i + Att , d^ A d^a - 

ye — e 1 ye — e 1 

ye*'' — e^^'i 

(C.23) 

Plugging the values of ^1 and ^2 in (IC.lip . and taking into account the proper normalization 

factors, we find exactly: 

--^sin^d0Ad0AW^ = 167r2(], (C.24) 



where Q is the one written in flC.4p with the following function S{ 



r : 



S{r) = Vl - e^'-^r, e(r - Vg) . (C.25) 



Notice that flC.25P is the limit of the function S{r) written in fl6.13p when r and r^ are large. 



This confirms our results of section |6l 

C.3 An example of a non-compatible embedding 

As we saw in the previous subsection, one has to work quite hard in order to obtain the 
smearing form Q from the microscopic average over a family of embeddings. Certainly, the 
trick described in section |5] gives a much faster and simpler way to get Q. One can wonder 
nevertheless about the reliability of the trick, since it assumes a given functional form of ^2, 
and the only unknown is the radial profile of the brane distribution S{r). 

In principle the trick can be run for any representative embedding. However, it is hard 
to think that any given family of embeddings, even if supersymmetric, generates (when we 
place fiavor branes along the embeddings of the family) a backreaction of the metric that is 
compatible with the initial ansatz we assumed for this metric, if this is not the most general 
possible. It seems nonetheless that the trick is able to detect this "compatibility property", 
and we present in what follows some arguments in favor of that. 

Recalling the discussion in section [5l the trick was to compute the effective radial action 
of the smeared brane distribution, and to compare it with Nf times the WZ effective ra- 
dial action of a single brane sitting at one of the embeddings of the family over which we 
smear. Both actions should be equal. The smeared action always contains two terms, one 
proportional to S{r), and another proportional to S'{r): 

^smeared ^ ^^^^) ^^^) ^ ^^^^) ^/^^) ^ ^^26) 
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where Fi and F2 depend on the functions of the ansatz. We conjecture that the way to 
detect if a family of embeddings generates a backreaction compatible with the ansatz is to 
take any representative embedding of this family and to compute its WZ effective radial 
action. We must then check whether or not the result depends functionally on the functions 
of the ansatz as in ( JC.261) . Let us assume that this is the case and that the effective WZ 
radial lagrangian density for the representative embedding is of the form: 

J^wt = Hr) G(r) + F,{r) H{r) , (C.27) 

where Fi and F2 are the same as in (1C.26I) and G{r) and H{r) are functions of r which 
do not depend on the functions of the ansatz. In order to verify that (1C.27I) is of the form 
flC.26|) one must check that: 

^ = H{r) . (C.28) 

dr 

If this were the case, we conjecture that the backreaction is compatible with the ansatz and, 

furthermore, that the profile function 5" is proportional to G. 

In section El we have worked out one example in which the compatibility condition was 
satisfied (see also appendix[D]). Moreover, in subsection IC.2l above. we have checked exphcitly 
with an independent calculation of Q that the trick gives the right result. In what follows, let 
us illustrate with an example the case in which the compatibility condition is not satisfied, 
and show with a microscopic calculation that indeed the resulting Q is incompatible with 
the ansatz for it. 

We choose to work again in the abelian background, since it is simpler and therefore the 
explanation will be cleaner. Let us focus on the following embedding: 

Ci = C^C4, C2=/x, (C.29) 

where the Q^s are the complex coordinates (1C.3I) and C and n are constants that we param- 
eterize as C = tan -^ e~'^'^° and /i = cos y e^'^°/'^ e*^. We can solve the embedding equations 
in ( 1C.29J) in terms of coordinates as: 



e = eo, = 00 , and e" cos I = e"', V' + = 2/3 . (C.30) 

It is easy then to compute the effective radial lagrangians of the smeared distribution and 
of a single brane extended along the embedding (IC.29|) . with the result: 



^smeared _ ^^ ^^ j.^^ ^2* ^ ^2k g ^ :_ g> 



e 



2g 



r^T = 27r Td5 e'* {e^' (l - e^"'-^') + Ae^^ e^'"-^') , (CSl^ 
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where we have assumed that Q should be as in (lC.4p . As we see, the e^^ term in the 
smeared lagrangian is not present in jC^f ^ (we have instead an e^'^ term) , so the compatibihty 
condition is not satisfied. 

Let us now check with a microscopic computation that, indeed, the family of embeddings 
generated by rotating (1C.29P with the SU{2)l x SU{2)r symmetry, generates an Vt that is 
not of the form of (1C.4I) . After using the relation ( IC.ll) . the family can be characterized by 
/i = 0, /2 = with: 

/i = aCi + K4 , 

f2 = K2 + lU + hix, (C.32) 

with |ap + |6p = 1 = |/p + \k\^. In this case, it is easy to perform the integral (1C.9I) by 
making use of the following two results: 

2? 
dzdz5^^Uwiz-W2) = —, — ^, (C.33) 

/ dx dy [x^ + y^ + aix + (3iy + 71) 5 [x^ + y"^ + a2X + I32y + 72) = 

= vr 7i - 72 H . (C.34) 



2 
The final result we get for the smearing form is: 

fi = —^ sin 6 d^^d<p^ ( (1 - e^^''-^'^) sin ^ d^ A d0 - 2e^''^-'^'dr A (d^ + cos 9 d<p)) , (C.35) 

and we see that this is clearly incompatible with ( ]C.4|) (the roles of {9, 0) and (^, 0) are 
exchanged in these two expressions of Vi). 

D Revisiting the massive KS background 

The fiavored version of the Klebanov-Strassler (KS) background |12], corresponding to a 
system of D3- and D7- branes on the deformed conifold, was obtained in ref. |T7] for the 
case of massless fiavors. This geometry was generalized in [23] to the case in which the 
fiavors are massive. In reference |23], in particular, a microscopic calculation of the D7-brane 
charge density was presented. In this appendix we will revisit this problem by applying the 
techniques developed in the present paper and we will reproduce and generalize the results 
found in [23]. 

The ten-dimensional metric in Einstein frame of fiavored KS mode is given by: 

ds^o = h'^ da;i,3 + h^ dsl , (D.l) 
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where h is a warp factor and dsg is the internal six-dimensional metric of the flavored 
deformed conifold [IT]. In order to write explicitly this metric, let us introduce the five 
one-forms g^ {i = i . . .5) as: 

, — sin 6 d(f) — cos ip sin 6d(f) + sin ip d6 o 

" = T2 • " 

— sin 6 d(f) + cos ijj sin 6d(f) — sin ip d6 4 



V-2 

g^ = d^ + cos 9 dcp + cos 9 1 
Then, dsl can be written as: 



d9- 


sini/j 


sin 6^ d0 — 


cos ip d9 






V2 




d9 + 


sinip 


sin ^ d0 + 


cos ^/^ d9 



V2 



(D.2) 



ds^ = r(r)e2G^M 



2\2 



(^^)^ + (^^) 



o2G2(r) 



(^ 



3^2 



(^ 



4\2 



jSGsCr) 

9 



(dr)2 + {g 



5n2 



(D.3) 
where the functions Gi{r) satisfy a system of first-order BPS equations and r(r) is the 
function: 

r(r) = ^^ . (D.4) 

cosh (r J + 1 

This background has a running dilaton $(r), a NSNS field ilf(3) = di?(2) and is endowed with 
RR one-, three- and five- forms F(i), F(3) and F(5). We refer to the article [23] for a complete 
account of this massive flavored solution. Here we will only need the expression of 5(2) and 
-F(i), which are given by: 



M 
B„, = - 



J(2) 



fg'Ag' + kg'^Ag' 



NfSjr) 
^(1) - 7Z — 9 



Att 



(D.5) 



where M is the fractional D3-brane Page charge, S{r) is the D7-brane charge density function 
(similar to the one used in the main text for the model dual to A/" = 1 SYM) and / and g 
are determined in terms of the dilaton $ as follows: 



/ 



r coth r 



2 sinh r 



-(coshr — 1) , 



k 



J, rcothr^ _ 

2sinhr ^ J k J 



Notice that dF(i) = —Q, where Q is the symmetry preserving D7-brane density distribution 
two-form, which is given by: 



n 



4:71 



Sir){g'Ag' -g'Ag') - S'ir)drAg' 



(D.7) 



The full background can be determined if the profile function S{r) is known (see section 
3.3 in [23]). In what follows we will show how one can find S{r) by applying the technique 
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developed in sections [MZl of the main text. In particular, we will be able to reproduce and 
simplify the results written in eqs. (2.14) and (2.15) of ref. |23] by using a holomorphic 
formulation similar to the one employed in the main text. 

The Kahler form of the transverse 6d manifold can be written as |23j : 

h-hj = e^'+'^'' (^g^ A g^ - g"" A g""^ - — dr A g^ . (D.8) 

The holomorphic variables Zi, z^-, ^3, z^ that we will use in this case are just the same as those 
defined in ( 14. Sp with the substitution r — )■ r/2. Similarly, the S'0(4)-invariant (l,l)-forms 
r/i, ri2 and 773 are defined in (14. lip . Their expressions in terms of r and the angles are just 
the ones written in (14.120 with the changes r — )■ r/2 and dr — )■ dr/2. One can easily show 
that it is possible to express J and B2 in terms of these forms as: 



J = -2i^-— r]i e^i+^2 + r/2 - — ^e 



« = -^^:^{'H^^^mU^i^m- (D.io) 



sinhr \ \sinh^r 9sinhr 

i^(2) = ,Me^ "^".^^"~% 3. (D.9) 

smh r 

The smearing form Q in this case, written in ( ]D.7p . can also be expressed in terms of the rji 
forms, namely: 

iNf S{r) / coshr \ iNf S'{r) 

27r sinh r V sinh^ r / 2tt sinh^ r 

The D7-branes wrap a non-compact four-cycle on the transverse space. For a supersymmetric 
configuration the corresponding eight-dimensional worldvolume is calibrated by the eight- 
form: 

/C = 1^ A [JAJ - e-*5(2)A5(2)]. (D.ll) 

It follows from (ID.9P that /C will also be expressed in terms of the S'O (4)-invariant (l,l)-forms 
■r]i , namely: 

/C = ^ Cij{r) d^x Ar],A rjj , (D.12) 

with Cij being the functions: 

2g2Gi+2G2 

sinh^ r 



4gGi+G2 / I 

C12 = ^7— cothre^^+^2 - -e^^^ 

sinh r \ 9 

9 G1+G2 / 1 

C22 = - "" . ( cothre«^+«^ - -e2«3 
sinh r 



(D.13) 



, e* (r cothr — 1) 



2 



C33 = M' — J— 

2h sinh r 



58 



The integration of the cahbration form over the D7-brane worldvolume will basically give 
the WZ term of the action. Our method to compute S{r) amounts to compare this term of 
the action for a smeared set of equivalent flavor branes with the same quantity evaluated for 
a single brane of the set. The former is given in terms of /C and Q as: 



^smeared = T7 / 6* /C A f] , (D.14) 

while the latter is just: 

S7t = T7 I ^*(e*/C) . (D.15) 

JMs 
Clearly, if the total number of flavor branes is Nf, one should have: 

5-^--^ = NfS'i^^t . (D.16) 

Let us show how (]D.16P can be used to determine S{r). In order to compute the WZ term 
of the smeared action (ID. 141) . let us define the functions Fi and F2 as follows: 



9 h 

Then, one can check by direct calculation that: 

e*/CAfi = ^e* \FiS + Fa^'ld^x A dr A d^G , (D.18) 

IGtt 

where d^G is the angular five-form: 

d^e = sin^sin^d^Ad0Ad^Ad0Ad?/^ . (D.19) 

Using these results in flD.14p . we obtain after integrating over the angular coordinates: 



^smjarcd _ -^j^^NfTr f d'o; dr 6* [F^S + F2S'] . (D.20) 

The next step is then to compute the pullback of the calibration form for one of the 
embeddings within the family used for the smearing in [23] . Let us consider the following 
particular representative embedding 



zi - Z2 = 2/i = 2/i e^" , (D.21) 

with fi being the modulus of fi, and a its phase. This embedding defines a four- dimensional 
submanifold inside the deformed conifold. Accordingly, we choose the complex coordinates 
Z2 and Z3 to parameterize it, where Z2 = Z2 + ft: 

z^ -1- u^ 
zi=Z2 + fi, Zi = -^ ^ . (D.22) 
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The relation between r and the complex coordinates follows from: 

2 coshr = y |z,|2 = 2|z2|' + 2^2 + \^^\2 ^ l^^-l-A'P _ ^^23) 

It is easy to compute the minimum of this expression. Let us denote by r^ the minimal value 
of r. Then, one has: 



2 cosh Tg = min s 2 1 ^2 1 + 2/i + min < | -23 1 + 






mm 



in{2|52p + 2/x2 + 2|z^-l-/i2|}_ = 2 [^^ + \1 + jl'^\) . (D.24) 



Varying the phase of fi, the minimum of this expression is always attained for a purely 
imaginary fi. Notice that when \fi\ < 1, it is possible to achieve Vg = 0. This matches the 
discussion in [23]. Moreover, for given /i, solving for the phase of jj, in (lD.24p . we get: 



_ / coshrg + 1 / l + 2/i^ -coshrg 
y cosh r, — 1 y 1 — 2//^ + cosh r^ ' 

which is exactly the expression (A. 10) of the paper [23] and tells us that actually the em- 
beddings flD.2ip are inequivalent for different values of the phase of ft. 

In order to compute the pullback of /C needed to calculate 5'^|*^, we have to calculate the 
pullbacks of the 4-forms rji A r]j. We find: 

I* iVi A r/j) = Aj (r) ^d52 A dl2 A dz^ A d^s , (D.26) 

where the only non-zero A- s are: 

All = 2 (coshr — yU^) , 

A12 = 2/i^ (cosh r + cos 2a) — sinh r , (D.27) 

^33 = — 2/i^ (1 + cosh^ r + 2 coshr cos 2a) . 
Now we have all the information to carry out our procedure. Let us write: 

i* (e* /C) = fe* V Cij{r)Aij{r)) -^ d^x A d52 A dia A d23 A dz^ . (D.28) 



Let us first compute the terms in parenthesis on the right-hand side of (]D.28p . From the 



values of the Cjj's written in (ID. 131) and those of the Aj/s displayed in ( ID. 271) one readily 
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gets: 



e" Y.^^,{r)A,,{ 



— e 



(sinh r — 2/i^ cosh r — 2/i^ cos 2a) Fi + 



2sinh'^r 

1 + cosh^ r + 2 cosh r cos 2a 



+ /i 



sinh r 



(D.29) 



where Fi and F2 are the functions defined in f lD.17p . This has a very suggestive form if we 
compare it with the resuh we obtained for the smeared action in ( 1D.20p . In order to continue 
with our calculation of S^^ ^ we use the fact that for any function F{r) we have: 



dz2 A dz2 A dzs A dz^ 



F{r) 



drJ{r)F{r) 



where the jacobian (which is computed below in subsection ID. 2^ is given by: 

J7'(r) = Svr^ sinhr . 
With these results we can write the action for a single embedding: 

sinh^ r — 2/i^ cosh r — 2/i^ cos 2a 



(D.30) 



(D.31) 



Qsinglc 



f 



-ATx'Tf / d^xdre 



sinh r 



F, + 



+ 2/i^ 



1 + cosh r + 2 cosh r cos 2a 
sinh^ r 



(D.32) 



By plugging ( ID. 201) and ( ID. 321) into ( ID. 161) . one immediately extracts the value of the profile 
function S{r): 

sinh^ r — 2jj? cosh r — 2/i^ cos 2a 



S{r) 



9(r — Tq 



sinh^ r 
An important consistency check of this identification is provided by the relation: 



d 
dr 



sinh^ r — 2/i^ cosh r — 2fi^ cos 2a 
sinh r 



2/i' 



1 + cosh r + 2 cosh r cos 2a 
sinh r 



(D.33) 



(D.34) 



Let us eliminate in (1D.32|) the phase a in favor of the minimal distance Vq. From ( ID.25P one 
can show that a is given by: 



2/i^ cos 2a = sinh^ r^ — 2yU^ cosh r^ . 

Then, one can rewrite the profile function in (ID.33P as: 

(coshr — coshr^) (coshr + coshr^ — 2/i^) 



Sir) 



sinh r 
61 



0(r - Vg) . 



(D.35) 



(D.36) 



Clearly S{rq) = as it should. Actually, near r = Tq the profile function behaves as: 

2(coshrq - 2/^2) 
S{r) -^- (r - Tq) , (D.37) 

which in particular means that S'{rq) ^ and, thus, S'{r) has a finite jump at r = r^ and 
the Ricci tensor will have the same type of singularity at the position tip of the branes. 

D.l Taming the threshold singularity 

As in section [7] we are going to find regular solutions by considering a brane distribution in 
which the position r^ of the tip of the branes is not fixed but distributed with a measure. 
First of all, let us define x and Xq as: 

X = coshr , Xq = coshr^ . (D.38) 

The density distribution function S{x) will be obtained by superposing the "elementary" 
distributions of eq. ( ID. 361) . One gets: 

Six) = r dxqpixq) i^-^.)i^^+^^.-^^^') ^ (D.39) 

Jxo ^ ^ 

where xq is the minimal value of x and p{xq) is the measure function, normalized as in (17. 3p . 
In the next two subsections we will consider two possible elections for the function p{xq). 

D.1.1 Phase average 



As mentioned above, for the embeddings (]D.21|) the position Xq of the tip of the branes is 



related to the phase a of the parameter fi by means of (ID. 240 . which now we rewrite as: 



Xq = p^ + V^l + /u4 + 2/^2 cos(2a) . (D.40) 

The prescription proposed in ref. [23] consists in averaging over all possible values of a with 
a flat measure p{xq)dxq ~ da. Actually, taking into account that Xq in (1D.40P depends on 
cos 2a, it is enough to consider a in the range < a < 7r/2. Let us denote by xq and xi the 
minimal and maximal values of Xq respectively. For /i > 1 (which we will assume in what 
follows) Xq and Xi are given by: 

xo = 2/1^ - 1 , xi = 2/i2 + 1 . (D.41) 

Then, the measure for this prescription is given by: 

2 r -I 

p{xq) dxq = — 0(xg — xo) — 0(xg — Xi) da , (D.42) 



TT 
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where the Xg's on the right-hand side should be regarded as the functions of a written in 

Let us compute the profile function 5* as a function of the radial variable x. From equation 
(ID.25p . we can see that the branes are lying in the range 2yU^ — 1 < Xg < 2yU^ + 1. We can 
distinguish three regions. 

Region I: x > 2/i^ + 1. In this region, the effect of all the branes (0 < a < 7r/2) can be 
felt. So, the range of integration in a will be from 7r/2 to 0. Therefore, we can write: 



S{x) = 3 /da " ~^~"^; /^ ^"^"" . (D.43) 



2 rt x2 - 1 - 2[i^x - 2/^2 cos 2a 
TT Jo x2 - 1 

The integration over a kills the term with cos 2a and one gets: 



X 



S[x) = 1 - 2/i2- , a; > 2/i2 + 1 (D.44) 

x"^ — 1 

which is exactly the expression (2.14) found in [23] for region I. 

Region II: 2/x^ — 1 < x < 2/i^ + 1. In this region, not all of the branes are contributing, 
only those which reach a minimum distance smaller than x. Since the minimum value of Xg 
is attained for a = 7r/2 it suffices with integrating in the region a{x) < a < ^, where a{x) 
is obtained solving in (ID. 251) . namely: 



Therefore, one has: 



, , (s + l)(l + 2/i2 -x) , , , 

a(x) = arctan J) '-) —^ '- . p.45 

^ ^ ' ^' a;- 1 l-2;u2 + a; ' ^ ^ 



V2 ^2_ 1 _o,,2. 



X - 1 - 2/i^a; - 2/2"^ cos 2a 
ua 

' a{x) 

Performing the integral one obtains: 



Six) = - da ^-^-^ . (D.46) 



n/ \ 2/7r /\\/, r. 2 ^ \ 2/i^ sin2a(a;) ,^ ^_, 

7r\2 /V x^ — 1/ TT x'^ — 1 

which can be rewritten as: 



S{x) =- 

TT 



2 I V (x-l)(l-2/i2 + x) n V x^-l' 



1 /(l + 2/i2_a;)(l-2/i2 + x) 

^2 V ^^^^1 
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(D.48) 



This looks much simpler than the expression written in [23]. However, one can check that 
both are equivalent. That is a consequence of the following non-trivial identity involving 
elliptic functions: 



where the function y42(x,/i^) has been defined in equation (2.16) of [23]. Remarkably, the 
charge density function S can be further simplified in this region, namely: 



S{x) = - 

n 



X \ , 1-X2 /(l + 2/i2-x)(l-2/i2 + x 



(^ ~ 2/i' ^y^ ) arccos [x + — -^) + 



x^ — \ J 2/^2 ^ V x^ — \ 



(D.50) 



Region III: x < 2^"^ — 1. Here no branes are contributing, so S{x) = 0. 

D.1.2 Rectangular measure 

We will now consider the possibility of distributing the branes homogeneously in Xq in some 
interval xq < Xq < xi, where xq = 2/^^ — 1 and xi — xq < 2. The corresponding measure 
is just the same as in (17. 4p with 6 = xi — xq and the resulting profile function S{x) can be 
straightforwardly obtained by performing the integral (]D.39|) . One gets: 



C/ \ 1 22; + Xq — 3/i 2 ^ ^ 

b[x) = - —— — — (X-Xo) , Xo<X<Xi, 

3 (xi — xo) (x"' — 1) 

(D.51) 

^, , 3x(x - 2/i^) + xqXi + (xo + xi) (3/i^ - Xo - xi) 

5(x) = ^^-^-^^ , x>x,. 

We can expand this result near x ~ xq, with the result: 

which leads to a regular background at x = Xq. 

In the particular case in which xq and x\ are taken as in the case of the phase average 
[i.e. given by (]D.41|) ). the above expression for S{x) reduces to: 

S{x) = '^l;^ ' f (x + 1 - 2/.^)^ , 2/^2-l<x<2/^2 + l, 

Q[x^ — 1 

(D.53) 

5(x) = 1-1^4^' x>2/.^ + l. 

3 x^ — 1 
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Now, near the endpoint x = Xq = 2/i^ — 1, the behavior of S{x) is simply: 

Si^) ~ A (^ - V + 1)' + • • • • (D-54) 



By plotting ( ID. 531) and ( ID.SOP together one concludes that the profile given by the simpler 



expression (ID. 530 is indeed very close to the one given by the more involved expression 



D.2 Details of the calculation of the Jacobian 



In this section we will detail the calculation leading to equations (ID.30P and flD.311) . Let us 



start by reparameterizing the complex variables Z2 and z^, which take values in the entire 
complex plane, in terms of their modulus and phase: 



Z2 = V^e'^' , Z3 = y^e'^^ . (D.55) 

In terms of those new variables, one can rewrite x = coshr, from flD.23p . as: 

1 + /i^ + 2/^2^3 + {u2 + UsY + 2/^2 cos(2a) - 2u2 (/i^ cos(2a - 2/32) + cos(2/32)) 

^ = o ' 

2m3 

(D.56) 
and the integral is now: 

X = / -. — r-dzs A dz2 A d^3 A dzs Fir) = / — dwa A dM3 A d/^s A d^s Fir) . (D.57) 

J Fsl J «3 

Now we replace u^ by x in the integral ( ID. 571) . The relation between these two variables is: 



M3 = a; - /i^ - M2 ± M , (D.58) 

where M is defined as: 

M = v^x^ - 1 + 2m2(cos(2/32) - x) + 2fi^{u2 -x- cos(2a) + U2 cos(2a - 2/32)) • (D-59) 



There are two different zones in the integral over ^3 in f lD.57|) . one for < M3 < ^3^0 



(where one uses the relation with the minus sign) and the other one for ^39 < "^3 (where one 
chooses the plus sign). M3 is the value of M3 for which x is minimum, which is given by: 



M3 = ^i + f^i + ul + 2fi^ cos(2a) - 2^2 {fi^ cos(2a - 2/32) + cos(2/32)) . (D.60) 

The value of x at that point is just ^3,0 = /i^ + ^2 + ms.o- Therefore, the integral over U3 can 
be written as: 



Us Jo Us i„3^ Us ^3^ M 
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Since the integrand in (]D.57p does not depend on the phase /Ss, we can write: 

l-2iT r-oo poo 1 

X = 47r / d/32 / du2 / —F{x) . (D.62) 

In order to perform next the integral over U2 we must take into account that x^^q depends 
on U2- By looking at (1D.25P one has that min{x3^o}M2 = ^g- Moreover x = x^^ for u = U2fi, 
with U2fl being: 



_ 1 1 + 2/^^x -x^ + 2/^^ cos(2a) 
U2 - U2,o - 2 ^2 _ 3, + ^2 cos(2a - 2/32) + cos(2/32) ■ ^ ' 

Therefore, one can exchange the order of the integrations over U2 and x as follows: 

dM2 / dx= dx du2 , (DM) 

J X'ifi J Xq Jo 

and the integral over U2 can be performed with the result: 

"^•« dM2 _ ^^2 _ I _ 2^23, _ 2/i2 cos(2a) 

M" ~ X - ;u2 - ^2 sin(2a) sin(2/32) - (/U^ cos(2a) + 1) cos(2/32) ' 

The only remaining integral to perform is the one on ^2- By employing the standard in- 
tegration techniques in the complex plane, this integral can be easily computed with the 
result: 



'" ^n A7T^x^-l-2f,^x-2f,^cos{2a) ^ 

,0 a; - /i2 - ^2 sin(2a) sin(2/32) - (/x^ cos(2a) + 1) cos(2/32) ^ ' ^ ' ' 

Remarkably, the right-hand side of flD.66P does not depend on x, a or /i. It follows that: 



/•OO /'OO 

I = 8tt^ / dxF{x) = 8tt^ / dr sinh rF(r) , (D.67) 

J Xq Jrq 

in agreement with (]D.30|) and ( ID.31|) . 
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